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I. 


INTRODUCTION 


This report covers the first year's work on the NASA/LERC 
research program on electrically pumped continuous wave (cw) excimer 
lasers. An excimer laser transition is a transition between two elec- 
tronic states of a diatomic molecule, the upper state being a bound 
state and the lower state being an unbound state. Such a transition 
will exhibit several advantageous features when compared with other 
types of lasers. These include 

« Broad band tunability (AX. ~1000 to 2000 A) 

e High cw saturation power (-50 kW/cm ) 

• Large high power system capability (low gain, 

fewer parasitic problems). 

The excimer laser transition is uniquely suited to meeting NASA's 
ultimate cw high power laser mission requirements. 

The primary objectives of the NASA program covered by this 
report are (1) to identify potential excimer molecules for visible and 
near visible cw laser operation based on electrical excitation, (2) to 
predict the specific operating conditions required for such laser action, 
and (3) to demonstrate cw laser gain on such a system m a laboratory 
experiment. 

To accomplish these objectives, the program was divided into 
four major tasks: 

Task I — Preliminary evaluation 

Task II — Analytic evaluation 

Task III — Laboratory apparatus 

Task IV — Molecular system evaluation. 

Specific objectives, procedures, and results for these four 
tasks will be presented m Sections II through V. Section VI presents 
preliminary conclusions and the work planned for the next year's 
follow-on phase of the program. 
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II 


PRELIMINARY EVALUATION 


The objective of this task is to identify at least three molecular 
excimer systems suitable for visible or near visible cw laser operation 
when excited by either a self-sustained or an electron beam- sustained 
discharge . 

Our procedure in addressing this task was to first select the 
three candidate molecules and to then select a cw electrical pumping 
technique appropriate for an experimental evaluation of two of the 
excimer systems under cw conditions. 


A. Selection of Candidate Excimer Molecules 

Consider again the basic definition of an excimer laser transi- 
tion as a transition that occurs between two electronic states of a dia- 
tomic molecule of which the lower state is a dissociative or repulsive 
state as in Figure II-l(a), In our selection of candidate cw excimer 
systems we broaden the definition of an excimer system to include any 
molecular electronic transition that exhibits a continuum red shift rela- 
tive to the transition between the associated parent atomic states, 
since it is this feature that is fundamentally responsible for the ability 
to maintain a cw population inversion on a broadened molecular elec- 
tronic transition. To see this, one assumes that the molecular states 
are populated m thermal equilibrium with respect to their parent atomic 
states m accordance with the gas temperature, while the population of 
the atomic states are in thermal equilibrium at the electron tempera- 
ture. The population ratio between the upper and lower molecular 
levels then can be written 


N 

u 



exp 


A(hv) 

T 

g 



( 1 ) 


where A(hv) is the red shift, hv^ is the atomic transition energy, and 

T and T is the gas temperature and electron temperature, respect- 
g e 

ively; we have omitted degeneracy factors for simplicity. Therefore, 
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an inversion will maintain, under this two -temperature equilibrium 
condition, if the ratio of gas to electron temperature is less than the 
fractional red shift, i. e. , 


T 

_£ 

T 

e 


Ahv 

hv 

o 


( 2 ) 


regardless of the bound or unbound nature of the molecular potential 
energy curves. An example situation, which exhibits the required red 
shift but does not involve a repulsive lower molecular level, is the 
dimer transition, illustrated in Figure II- 1(b), where the red shift is a 
result of the shift in the mternuclear separations of the potential mini- 
ma of two bound molecular states. 


4053-1 4053-2 




o b 

Figure II- 1. Excimer and dimer potential energy curves, (a) Excimer. 
(b) Dimer. 
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The two -temperature equilibrium condition is a sufficient but 
not a necessary condition for achieving laser gam on red shifted exci- 
mer and dimer transitions. For example, excimer and dimer laser 
action demonstrated on the vacuum ultraviolet rare-gas excimer and 
the dimer transition, respectively, does not occur under two- 
temperature equilibrium conditions. It is still the red shift, however, 
that is primarily associated with the ability to achieve a population 
inversion in these cases. 

The search for potential excimer molecules began with the 
observation that every element in the periodic table will probably form 
at least one excimer or excimer-like (m the sense defined above) laser 
transition either m a homonuclear or heteronuclear diatomic molecule. 
Thus, our preliminary evaluation of potential excimer molecules first 
involved various processes of elimination applied to the elements of 
the periodic table (Figure II-2). 

The first criterion applied in the elimination procedure is 
somewhat arbitrary but motivated by the desire for expediency. Since 
electrical excitation is to be used, it is required tbat the active medium 
be m the gas phase. In particular, we consider only those elements 
which have a vapor pressure greater than one Torr at temperatures 
below 500°C, and are reasonably abundant and nonradioactive . This 
leaves us with the following elements (see Figure II - 2 ) : 


Group 1 
Group II 

Transition Elements 

Group III 

Group IV 

Group V 

Group VI 

Group VII 

Group VIII 


H, Na, K, Rb, Cs 
None 

Zn, Gd, Hg 

None 

None 

N, P, As 

O, S, Se 

F, Cl, Br, I 
He, Ne, A, Kr, Xe 


A further narrowing down of the group of candidate excimer sys- 
tems was then achieved through qualitative considerations on the nature 
of the known excimer or excimer-Like transitions on both the heternuclear 
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4388-5 



Figure II-2. Evaluation procedure applied to the periodic table 






diatomic molecule formed by the association of the element with a rare 
gas atom as well as on the homonuclear diatomic molecule. Figure II-3 
shows potential energy curves of a heteronuclear (excimer) and homo- 
nuclear (dimer) candidate molecule representative of each group, 
together with a brief listing of problems and advantages associated with 
the use of the system as a cw excimer laser. References consulted for 
the survey are listed in Table II- 1. A brief narrative summary of our 
considerations of these systems follows. 

4388-6 

EXCIMER (HXe) 


Problems: Wrong 

Wave length 

Candidate Systems No 


DIMER (H z ) 

Problems: 

No net gam due to photoionization 
of upper laser level 

Candidate System ?: No 


interimuclear separation do -2 cm* Figure II- 3 . Group 1 Hydrogen. 
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Table II- 1. References Consulted for the 

Preliminary Evaluation Task 


EXIiVIE R CANDIDATE 

REFERENCE 

HXe 

1,2 

CM 

X 

3 

Kx e 

4 

K 2 

5 

HgXe 

1 

Hg 2 

6,7 

Xe 0 

9 

Xe F 

10 

Xe 2 

11 


T 1 937 


Group I; Hydrogen — The bound-free excimer transition on 

the hydrogen-rare gas molecule, which correlates with the resonance 

transition, will lie in the vacuum ultraviolet and therefore is not a 

o 

candidate system. The well-known bound-free transition at 3500 A of 

H_ has been investigated as a possible laser transition by Palmer 
^ 3 

several years ago. It was found in that investigation that photoioniza- 
tion transitions from the upper laser level will not permit net gain to 
be realized on this system. The associated bound-bound transitions of 
are in the vacuum ultraviolet and involve well depths far too deep 
for the transition to be broadened into a continuum as is necessary to 
qualify it as an "excimer-like" transition. Hydrogen was therefore 
judged as not a good candidate element for this program. 

Group I: Alkali Metals — The alkali metals are known to be 
excellent candidate excimer laser systems m both the A-X bound- 
unbound transition on the alkali-rare gas molecules and the A-X bound- 
bound transition of the alkali dimer molecules (Figure II-4). The 
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4388-7 



8 13 18 


Advantages; 

Good wavelength (near IR — 
visible) 

Strong oscillator strength 

No competing photoexcitation 
processes 

Low pumping threshold (flash- 
lamps, discharges) 

Problems: 

Weak upper level binding 
(-0.1 eV) 

Thermal bottle-necking and 
dimer absorption at moderate 
temperatures 

High pressure alkali-rare gas 
mixture, difficult to handle 

Candidate System: Yes 



8 13 18 23 28 


Advantages: 

Same as excimer except: 
potentially broader wave- 
length range and operation 
at lower pressure 

Problems: 

Same as excimer except: 
can operate at slightly higher 
temperature s 

Candidate System: Yes 


Figure II-4, Group 1 — Alkali metals: Na, K, Rb, Cs. 
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dimer transition is red shifted as required and is broadened into a con- 
tinuum at buffer gas pressures greater than or near one atmosphere. 
The excimer laser features of these systems are discussed in detail 
m Refs. 5 and 12, Possible problems associated with the use of these 
systems for cw operation are thermal dissociation of the upper levels 
resulting from the relatively shallow well depths and discharge insta- 
bilities caused by the metal vapor gradients. The K-Xe/K^ system 
was selected as the first of the three candidate cw excimer systems to 
be chosen from the preliminary evaluation. Potassium was chosen in 
favor of the other alkalis because it requires the lowest oven vaporiza- 
tion temperatures of the alkalies whose excimer and dimer transitions 
both lie within the range of a silicon detector. Xenon is chosen in 
favor of the other rare gases because it forms the deepest excimer 
state well depth (resulting in the farthest red shift). 

Zn, Cd, Hg Transition Element Subgroup — These elements 
have bound free transitions to the ground state on both the heternuclear 
rare gas molecules and on the homonuclear dimer molecules (Figure 
II- 5). 

In particular, Hg_ has both a visible excimer transition at 
o z o 

4850 A and an ultraviolet excimer transition at 3350 A which have been 

£ 

recognized for some time as potential laser transitions. Gain mea- 

o 

surements carried out on the 4850 A transition have revealed the same 

problem with this transition as with the excimer transition, i. e. , 

no net gain due to vapor state absorption. ^ On the other hand, the 
o 

3350 A transition has now been shown to exhibit net gain both with 
probe laser gain measurements as well as with a recent demonstration 

* 3 

of oscillation under optical pumping by a laser. The Hg„ 1 — y* 

° °2 u ^g 

33 50 A excimer transition was identified as another candidate cw 
excimer system. 

The corresponding heternuclear Hg-rare gas excimer transi- 
tion is not well studied. One can anticipate a weak oscillator strength 
for the transition and a weak excimer state binding a relatively short 

V. Schlie, private communication. 
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transition wavelength, probably about 27 00 A, which will elevate the 
pumping power required for threshold. Thus the heteronuclear Hg-rare 
gas molecule is not admitted as a candidate system. 


Group V and VI — The excimer and dimer transition involving 
the Group V elements are similar to those involving the Group VI ele- 
ments (Figure II-6). Of the elements m these groups, oxygen has 
recently been particularly well studied as an excimer laser. The 
recently demonstrated XeO e-beam pumped excimer laser operates on 
a bound-bound molecular transition which correlates with the atomic 
oxygen auroral transition. While this system may be well suited to 
large energy storage because of the long lifetime, it is poorly suited to 
cw application because of the evident self-terminating nature of the 
transition. The transition correlating with the resonance transition on 
both the heternuclear excimer and homonuclear dimer molecules are 
in the vacuum ultraviolet and are also unsuitable for the present program. 


EXCIMER (HgXe) 



4388— B 

DIMER (Hg 2 ) 



Problems: 

Weak oscillator strength 
Short wavelength 
Weak upper state binding 
Lack of data and theory on 
potential curves 

Candidate System: No 


Advantages: ^ 1 

Measured gain on 1 — 1 ► S 

° u g 

transition & 

Good wavelength 
Problems: 

Weak upper state binding 
Candidate System: Yes 


Figure II- 5. Transition element subgroup: Zn, Cd, Hg. 
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EXCIMER (X e 0) 



Problems: 

Extremely weak upper state 
binding ~0. 05 eV 

Evidently a self terminating 
transition 

Colli sionally induced dipole 
moment requires high pump- 
ing and high pressure 

Candidate Systems: No 


DIMER (0 2 ) 



1 2 3 


r 10~®cm 

Problems: 

Lower level vibration structure 
Too large: cannot broaden 

Transition into a continuum 
and collisionally deactivate 
lower laser level 

Candidate System: No 


Figure II-6. Groups V and VI; N, P, As, O, S, Se. 

Group VII — The halogen-rare gas excimer molecules 

appeared to be one of the most promising systems of all for use as a 

cw excimer laser (Figure II-7). Strong ionic-bonded excimer levels 

are formed by the association of a metastable rare-gas atom with a 

ground state halogen and the strongly allowed transition to the repuls- 

o 

ive or weakly bound ground state is near 3000 A, The deep excimer 
state well (~4 eV) will permit the use of high cw input power densities 
without upper state dissociation as occurs with the alkali systems. 

Also, because F, Br, and Cl are room temperature gases, the 
experimental apparatus is enormously simplified. The XeF excimer 
transition at 3450 X was chosen as the candidate excimer molecule 
from this class because it has the longest wavelength. The halogen 
dimer bound- bound transition correlating with the resonance transition 
is not suitable because of the short wavelength and deep well depths, 
with vibrational level spacing too large to be broadened into a continuum. 
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XeF 



Advantages: 

Ideal wavelength 
Strong oscillator strength 
Strong upper state binding 4 eV 
Fast excimer association rates 
10" ^ cm^ sec" 

Discharge pumping 
F, Be, Cl gasses at STP 

Problems: 

Discharge stability*? 

Candidate System: Yes 

Figure II-7. 



Problem: 

Same as for Group V and VI 
dimers 

Candidate System: No 


Group VII: F, Bi, Cl, I. 


Group VIII — The rare gas dimer bound-unbound transitions 
are well-known laser transitions, but they lie in the vacuum ultra- 
violet thus requiring very large pumping thresholds, and thus are not 
candidate systems for this program (Figure II-8). 

Having identified four specific candidate excimer molecules for 
detailed theoretical modeling our next task was to assemble prelimi- 
nary quantitative data on the laser transition. The four candidate exci- 
mer transitions selected and their approximate wavelength are 


(1) 

KXe • It 1 / Z) 3/2 S l/2 

8500 1 

(2) 

K 2 ^u* — 

10400 A 

(3) 

Hg 2 A 3 Iu 

3350 A 

(4) 

XeF S l/2 S l/2 

3450 % 
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A compilation of pertinent quantitative data, including a set of potential 
energy curves for these transitions, are presented in Figure II- 9 
through II- 11. References used for these data are the same as those 
listed in Table II- 1. 


B. Selection of a Cff Pumping Technique 

In selecting a cw pumping technique for excimer lasers it is 
important that one appreciates some rather general constraints con- 
cerning cw operation of an excimer laser. 

The stimulated emission cross sections for visible-near visible 

- 17 2 

excimer lasers average on the order of 3 x 10 cm . To sustain a 

practical laser gain of 1% per centimeter required for laser action it 

.14 

is necessary to sustain an excimer level population density of ~3 x 10 
- 3 

cm . Visible transition excimer levels have a radiative decay rate of 
~3 x 10 7 sec . Discharge pumping efficiencies, from recent experi- 
ments, seem to be running near ~1%. Thus, assuming a 1 to 3 eV 

4 5/3 

photon, a minimum cw input power density of ~10 to 10 W/cm must 
be provided to the excimer system for it to oscillate near threshold. 


4388-11 


Xe + Xe 




Xe(L) + Xe 


Problems: 

Wrong Wavelength 

Candidate Systems: 
No 


INTER NUCLEAR SEPARATION, R 


Figure II-8. Group VIII: He, Ne, Ar, Kr, Xe. 
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ENERGY x 10 


KXe 


4388-13 




EXCIMER FORMATION REACTION 

K(P 2^3 1 / 2 ^ ^ ^ KXe{^7T^2 3 / 2 ) Xe 

{RATE = 8 2 . 10 - 32 cm 6 sec -1 ) 

LASER TRANSITION 2rr y2 3!2~~ 2 ^V2 
TRANSITION WAVELENGTH X = 8200 - 9000 A 
UPPER LEVEL LIFETIME T ~ 2 7 • 10 -8 sec 


STIMULATED CROSS SECTION 



= 4.10 -16 cm 2 


V ' V o 

FRACTIONAL RED SHIFT ~0 092 


{EQUALS MAXIMUM Tg/Te FOR 
2-TEMPERATURE INVERSION) 


K{P 1/2 ) + K{S 1/2 ) + X — K^ 1 Ij+X 
(RATE = 8 • 10 -30 cm 6 sec -1 ) 

1 y ^ ^ 1 + 

X ~ 9300 -10,500 A 
T ~ 2 5 • 10 -3 sec 
cr=sl • 10 -16 cm 2 


V - V 


o 


0 19 


Figure II- 9. Laser Transition Data: Kxe, K^. 


21 





DIMER (Hg 2 ) 


4388-14 



Excimer formation reaction: 

Hg( 3 P x ) + Hg{ X S o ) + X-Hg 2 (A 3 l u ) + X 

(Rate * 1 x 10 30 cm^ sec" 1 ) 
Laser transition: 

Transition wavelength: 

\ « 3000 - 3600 1 

Upper Level Lifetime: r = 0.22 usee 

U n 

Stimulated emission cross section; 

I O O 

<r= 0. 97 x 10" cm 
Fractional red shift: v v o - q 34 


Figure II- 10. Laser Transition data: Hg_. 


XeF 


4388-15R1 



EXCIMER FORMATION REACTION 

Xe( 3 P 2 ) + F 2 — XeF( 2 S 1/2 ) + F 
(RATE = 2 -10 -9 cm 3 sec -1 ) 

LASER TRANSITION XeF( 2 2 1/2 ) — XeF( 2 2 1/2 ) 

TRANSITION WAVELENGTH X= 3250 -3500 A 

UPPER LEVEL LIFETIME r u = 50 ■ 10' 9 sec 

STIMULATED EMISSION CROSS SECTION 

cr = 05 • 10~ 17 cm 2 


v-v Q 

FRACTIONAL RED SHIFT = 0 58 


Figure II— 11. Laser transition data: XeF. 
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If the excimer laser is to be run in a truly steady- state cw 
operating mode, the power density (less the —1% laser power extracted) 
must also be continuously removed from the system and the gas 
temperature of the active medium must be kept below ~2000°C if the 
population inversion is to be sustained. 

The only three basic mechanisms for removing heat from the 
system are (1) radiation, (2) convection, and (3) conduction. We con- 
sider now the requirements for each of these mechanisms to be applied 
to a cw excimer laser. 

1 . Radiation 

For radiative cooling to reach the minimum required 
4 2 

rate of 10 W/cm (assuming a transverse dimension of ~0. 1 to 1 cm), 
the active medium would have to have an effective blackbody tempera- 
ture of ~1 eV. An excimer system m thermal equilibrium at this 
temperature will not exhibit a population inversion. 

2. Convection 

Convective cooling using gas flow at or near atmospheric 

pressure is commonly used to remove heat from high average power 

CO., lasers such as the Hughes Peacemaker laser. In these devices 

the temperature is controlled simply by the residence time of a volume 

element of gas together with its specific heat: (3/2) NT = p • ( d / v jq ow )> 

where N is the gas density, p is the input power density, V £^ QW IS the 

flow velocity, and d is the discharge dimension along the flow direction. 

It is easy to see that this technique is impractical for excimer lasers 
5/3 

requiring a 10 W/cm power dissipation. With an input power of this 

19 

magnitude at atmospheric pressure, sonic fLow conditions (N = 3 x 10 
3 5 

cm" , v.. = 10 cm/sec), d would have to be Si mm to keep T < 

flow ' 

1000°C. If the lower end of the above threshold pumping range is 
relevant (~10^ W/cm\ then convection may be practical. (The 
alkali excimer systems are expected to have pumping requirements 
near this latter value. ) 
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3. 


Conduction 


Thermal conduction of the heat to a wall m a static gas 

4/2 

and subsequent removal by flowing coolants can achieve 10 W /cm . 

This will be obtained for a helium gas mixture (independent of pressure) 

across an assumed thermal gradient of ~1000°C, which is close to the 

upper limit defined by thermal bottlenecking of the population inversion 

5 . 3 

and by the integrity of the wall material. For the 10 W/cm pumping 
requirement, sur face-to-volume ratios of ~10 will suffice for this mode 
of cooling. In spite of the obvious nonscalability of the thermal conduc- 
tion mode, the simplicity and versatility of this mechanism motivated 
the choice of a capillary discharge configuration to be used for the first 
year phase of the program. Also, capillary discharge pumping of CO^ 
waveguide lasers is now routinely used under cw conditions at pressures 
comparable to the pressures required of an excimer laser although 
under significantly lower power loading. Capillary discharges are also 
used for argon ion lasers under comparable power loading but at pres- 
sures lower than required for an excimer laser. 

While full cw operation, i.e. for durations of minutes, is not 
intended, experiments with pulse lengths up to a millisecond can be 
carried without the need for flowing wall coolants, and this will enable 
one to demonstrate cw operation conditions as far as the excimer laser 
kinetics is concerned. 
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III. 


THEORETICAL EVALUATION 


This section summarizes modeling procedures and results for 
discharge pumping of the XeF, KXe, and excimer systems. In the 
case of XeF, a comparison between the results of modeling and experi- 
ment is also presented. 


A. Modeling Procedures 

In each excimer system the model follows, by means of numeri- 
cal integration of the appropriate rate equations, the time evolution m 
the discharge of several excited and ionized species, the electron and 
gas temperatures, and the net small signal gam coefficient near the 
peak of the excimer gam band. The time dependence of the circuit 
currents and voltages is also included with a series resistance, a 
series inductance, and a storage capacitor as the modeled circuit 
elements. 

The average electron energy is computed by solving the electron 
energy equation. 


J ' E disch 



inelastic collision rate (Te))^ x E^ 


1 


+ (elastic collsion rate) x 2 



x Te 


( 3 ) 


where J is the current density, E .. , is the discharge electric field, 

1 disch^ 

E. is the ineLastic energy loss for the i process, and M, , and M are 
the masses of the major component species and electron, respectively. 
At present, the resonance or metastable level pumping is assumed to 
be the dominant inelastic process contributing to the energy equation. 
For this reason, together with the fact that the molecular association 
rates which pump the upper laser levels dominate the loss of the reso- 
nance level population, the predicted efficiencies for converting dis- 
charge power to laser power is essentially the quantum efficiency of 
the transition that is lasing. A Maxwellian energy distribution is 


25 



assumed. The Maxwell distribution should be a fair approximation 
because of the dominance of two-step ionization processes over single 
step ionization under the considered operating conditions. Only the 
latter process requires electrons with energies in the tail of the dis- 
tribution where a large reduction of the number of electrons from a 
Maxwellian distribution is known to occur. 

The electron collisional rate constants are expressed in a form 

that results from assuming a linear rise of the cross section with 

1 3 

energy above threshold : 

(sT Te ) 3/2 -Pl-o/Te) (l + ^gr), 

£ o (4) 

where e !s the threshold energy for the excitation or ionization 
process, dtr/dej^ is the rate of rise of cross section with energy 
above threshold, ^nd m^ and e are electron mass and change. 

Also, in solving eq. (3) we use an expression for the drift 
velocity obtained from Ref. 13: 

V drift = < 2 T '/ M M> 1/2 • < 5 > 

The gas temperature is calculated by assuming that all of the 
discharge power goes into gas heating and that heat is removed only by 
thermal conduction to a wall (assumed cylindrical). 


R(Te) = 


6 x 10 dtr 
(e/m)/ it de 


dT G 

dt 


J ' ^7(3/2 [M] ) - (thermal conductivity) 


4.8 

diameter 


X < T G " T wall> * 


The discharge circuit response is governed by the familiar circuit 
equations . 
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R x I 


(7) 


V dxsch = v c (t) ' L X f - 


V c<‘> = V charg e -VC flit , 


( 8 ) 


where ^disch is the discharge voltage, V is the capacitor voltage, 
^diar^e 13 * nit *- a ^ charging voltage, I is the discharge current, and 
R,L, and C are the series resistance, inductance, and storage capaci- 
tance, respectively. 


1 . XeF System 

The plasma, gas kinetic, and radiative processes included 
in the rate equations governing the species concentrations are presented 
in Tables III- 1 through III-3, together with source references and m 
the computer code listing presented m Section 1 of the appendix. (In 
the listing, a B at the end of the label refers to the back reaction which 
is computed by detailed balancing using the appropriate equilibrium con- 
stant or Boltzmann factor for the reaction. ) 

The small signal gam coefficient for the XeF system is computed 
as follows: 


gain = <r atim x ([XeF-] - [XeFg] x exp (- 0. 074/T g )) (9) 

where, 

(T (P) = ", X ( A ' X \ X A . , (10) 

Stlm 4x 2 c l AX (P)/ 6 “ 

and Ak is the measured half width of excimer band as a function of 

o 

pressure, k is the wavelength (3450 A for XeF*), A f s the A coeffi- 

cient for the transition, [XeF*] and [XeF ] are the upper and lower 
laser level concentrations, and the Boltzmann function accounts for gas 
heating effects on the population distribution in the lower level. 

The plasma, gas kinetic, and radiative processes modeled for 
KXe/K 2 excimer systems are listed in Tables III-4 through III-6, 
together with source references and in the computer code listing 
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Table III- 1 - Plasma Processes Modeled, for the XeF Excimer System 



PROCESS 



REFERENCE 


e + Xe Xe + e 


0 3 x 10 -10 cm 2 eV -1 
0 05x 10 _ ' 16 cm 2 eV~ 1 


Rough Estimate 


e + Xe Xe** + e ■ 


e + Xe Xe + + 2e 


0 11 x 10 -16 cm 2 eV -1 


e + Xe 4 -*Xe + + 2e“ 


e + Xe**-*Xe +2e' 


5 x 10“ 16 cm 2 eV~ 1 

—16 o _i 
5x10 cm^eV 1 


Rough Estimate 


Rough Estimate 


e + Xe 2 -*■ Xe* + Xe 


(a =) 0 346 x 1 0 — ° (T , v , )“ 


cm® sec ^ 


R4MX 


e + XeM -» Xe* + M 


(a =) 0 346x 10-°x (T e (eV) )' 
cm® sec~^ 


Rough Estimate 


Corresponding Processes for M species = Helium 


e +M2M 4 + e' 


0 025 x 10~ 17 cm 2 eV -1 


e +M^_M**+e' 


0 025 x 10 — 1 7 cm 2 eV -1 


Rough Estimate 


e + M -> M + + 2e 


0 08 x 10 -17 cm 2 eV -1 


e~ + M* ->■ M + 2e~ 


5x 10~ 16 cm 2 eV -1 


Rough Estimate 


e + M** -*• M + + 2e — 


e +M 2 + ->M* + M 


5x 10 -16 cm 2 eV -1 


(a =) 0 29 x 10 -8 {T e )" 1/2 
cm® sec - ^ 


Rough Estimate 


e + M -»-e + M 


{a =) 5 5x 10~ 16 cm 2 
{elastic collisions) 


e + F 2 -*• F + F 


()3=) 9 1 x 10" 10 (T e r 2 16 
exp {—0 165/T e ) cm® sec - ^ 


e + Xe F XeF* + e 


{R, =}2x10 8 cm® sec ^ 


e + Xe + wall -»• Xe 


(ambipolar diffusion) 

Da = 2x 10 3 x (3x 10 16 /[M] > 
v/ 0-03/Tg Te 
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Table III-2. Gas Kinetic Processes Modeled for the XeF Excimer System 


LABEL 

PROCESS 

RATE 

COEFFICIENT 

REFERENCE 

G3 

M + + M + M Z M 2 + + M 

3 7 x 10~ 31 xT g ’ /2 

All gas kinetic rates were obtained 

G3X 

Xe + +Xe + M^Xs 2 + + IVt 

12 x 10~ 31 xT g ' /2 

from Ref 17 or from scaling 
arguments applied to rates given 

G3MX 

Xe h + M + M Z XeM + + M 

5 8 x 10 — 31 xT g ’ /s 

in Ref 17 

G5 

M* 4- M 4- M Tm 2 + + M 

12 x 10~ 33 Tg 54 


G5X 

Xe* + Xe + WlTXe 2 ‘ + M 

15 x 10 _32 T g ' /2 


G5MX 

Xe' + M + M-^XeM' + M 

47x10~ 33 T g ' /2 


G13 

M * 4- Xe ~ Xe + + e — + M 

29 x 10- 10 T G * 


G14 

M 2 ' + Xe — Xe + + e — + M 

29 x 10 _10 T g ' /2 


G18X 

Xe" 4- F 2 XeF* 4- F 

34 7 x 10"" 10 xT g ' /2 


G20X 

Xe 2 * + Fj ■+ XeF* + F 4- Xe 

35 x 1C“ 10 xT g ' /2 


G22 

M + 4- F — -* M* 4- F 

5x 10~ 7 T g ' /j 


G22X 

Xe + + F -> Xe * + F 

17 x 10 -7 Tq 11 


G23 

M + +F"M-*MF‘ + M 

3 x 10 -25 T G y ‘ 


G23X 

Xe* + F + ft/1 ^ XeF * + M 

12 x 10“ 25 T G V ’ 


G24 

M 2 + 4-F~-»M 2 '4- F 

5 x 1 0 -7 T g ’ /2 


G24X 

Xe 2 4* F “*■ Xe 2 * 4- F 

17 x 10 — 7 T g ’ /2 


G25 

M 2 + 4- F - ^ Xe 2 * 4- F 

3 x 10" 25 T g ’ /2 


G25X 

Xe 2 + 4- F~ 4- M XeF* 

12 x 10 _25 T g ’ /2 



4- F4-M 

46 x 10“ 10 Tq' 7, 


G27X 

XeF* 4- F 2 ~ XeF + F 4- F 


G31 

M* 4- M* -» M 2 + 4- e — 

10 x 10 _10 T g ’ /2 


G31 X 

Xe‘4-Xe‘-^Xe 2 f + e - 

30 x 10'' 1 °T g ' /2 


G31MX 

Xe* 4- M* -> XeM + 4- 



G32 

M*4-M 2 *-^M 2 + 4-M4-e~ 

10 x 10 _10 T g ' /2 


G32X 

Xe* + X &2 Xe 2 * + Xe + e 

20 x io _10 t g 1/2 


G33 

M 2 * 4- M 2 * —■ M 2 + 4- 2M 4- e “ 

10 x 10~ 10 T g ,/4 


G33X 

Xe 2 * + Xe 2 * -* Xa 2 + 4- 2Xe 

30 x 10 _10 T g ’ /2 
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Table III- 2. Continued 


LABEL 

PROCESS 

RATE 

COEFFICIENT 

REFERENCE 

G33MX 

XeM* + XeM’ -*■ XeM + + M +Xe + e~ 

20 x 1 0 -1 0 T g ’ /j 


G34 

M* + XeF* -*■ XeM* + F + e — 

10 x 10“ 10 T g % 


G34X 

Xe* + XeF* -*■ Xr 2 + + F + e~ 

30 x 10 -,0 T g I/4 


G35X 

Xe 2 * + XeF* -*► Xe 2 + + XeF + e~ 

30 x 10 _,0 T g ’ /2 


G35MX 

XeM* + XeF* ->• XeM + + XeF + e" 

20 x 10 _,0 T g ' /2 


G40 

f + f + m-»f 2 + m 

1 x 10~ 32 T g ' /j 


G41 

F + Xe + M^FXe + M 

2 x 10~ 32 T g ' /j 


G44 

F + Xe* HI-* XeF* + M 

6 x 10 _32 T g ' /2 



Equilibrium Constants for Processes 



3X, 5X,and 41 


K3X 

1 x 10 — 23 exp (0 6/T q ) 

Rough Estimate 

K5X 

1 x 10 — 23 exp (1 5/T g ) 

Rough Estimate 

K41 

1 x 10 — 23 exp (0 074 /Tq) 

Rough Estimate 
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Table III-3. Radiative Processes Modeled for the XeF Excimer System 


LABEL 

PROCESS 

RATE 

REFERENCE 

AX2 

Xe**^Xe +hu 

4 x 1 0 7 sec - ^ 

Rough Estimate, 

AXED1 

Xe 2 * -» Xe + Xe + hu 

4 x 10 7 sec - ^ 

17 

AXEF 

XeF* XeF + hu 

2 x 10 7 sec - 1 

17 

AMX 

XeM* Xe + M + hu 

10 x 10 7 sec - ^ 

Rough Estimate 

am 2 

M** -?■ M + hu 

32 x 10 7 sec -1 

Rough Estimate 

AMD1 

M 2 * JV1 + M + ho 

32 x 1 0 7 sec -1 

Scaled from Xe 2 * 
value in Ref 17 


T1945 


Table III-4. Plasma Processes Modeled for the KXe/K, Excimer System 




do 


LABEL 

PROCESS 

dt e0 

REFERENCE 

RPI 

e - + K -s- K + + 2e — 

0 5 x 10 -16 cm 2 eV -1 

14 

RPI 

e - +K^K* + e - 

6 x 10 -15 cm 2 eV -1 

14 

RPII 

e~ + K* - K + + 2e — 

5 x 10 -16 cm 2 eV - ' 1 

Rough Estimate 

RECP 

e + K + + e K* + e 

(a =) 2 28 x 10 -26 (T g ) -4 39 

fi 1 

cm sec 

15 (cesium) 

RECPX 

e - + K Xe + ->• K* + Xe 

(a. =) 3 3 x 1 0 — 5 (T e ) -0 67 

3 -1 

cm sec 

Rough Estimate 

RECPD 

e — + K 2 + -9- K*+ K 

(a =)33x10 -5 (T e ) — 9 67 

3 -1 

cm sec 

Rough Estimate 

CEL 

e — + Xe - 9 - e — + Xe 

( 0 =) 55x10 -16 cm 2 

14 
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Table III-5. Gas Kinetic Processes Modeled for the 
KXe/K^ Excimer System 


LABEL 

PROCESS 

EQUILIBRIUM CONSTANTS 
OR 

RATE COEFFICIENT 

REFERENCE 

KXEX 

K + 2 Xe ^ K Xe + Xe 

2 2 x 10“ 23 exp {-0 05/T G )cm 3 

Rough Estimate 

KXD1 

K + K + Xe Z K 2 + Xe 

1 8 x 10~ 22 exp (0.56/T G )cm 6 

5 

KAXI 

K* + 2 Xe +1 K Xe* + Xe 

2 2 x 10“ 23 exp (0 074/T G )cm 3 

18 

GPD 

K* + K + Xe Z K/ + Xe 

8x10 cm sec 

5 

KAXI 

K + + 2 Xe Z K Xe* + Xe 

2 2 x 10~ 23 exp (0.074/T g ) cm 3 

18 

KADI 

K + + K + Xe Z K 2 + + Xe 

6 5 x 10 -23 exp (0 735 /Tq) cm 3 

Rough Estimate 

CD IF 

Heavy particle elastic 
cross section 

1 x 10~ 15 cm 2 

Rough Estimate 
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Table III- 6. Radiative Process Modeled for the KXe/K^ Excimer System 


LABEL 

PROCESS 

RATE 

REFERENCE 

AT 

K*^K + hv 

Trapped Rate 

20 



(computed) 


AP 

KXe* -> KXe + ho 

3 69 x 1 0 7 sec --1 

19 

AP 

K 2* ^ K 2 + ho 

3 69 x 1 0 2 sec -1 

Rough Estimate 


T1943 
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presented in Section 2 of the Appendix, In this system, under the 
conditions we consider ([Xe] ' 1 x 10 cm ) the gas kinetic rates 
dominate over electron collisional rates in controlling the population 
of the molecular states. We therefore assume thermal equilibrium at 
the gas temperature to be maintained between the molecular states and 
their dissociation products. An exception is the A state of K for 
which radiative losses are included since they can compete with the 
dissociation rate. 

The small signal gain coefficient for the KXe/K system is 
computed from the quasi-state theory of line broadening as in Ref. 21: 


exp (-V A ^ (X)/t ) 

«<M |[S*-] * (K ° 

' eq A, Ex 


1/T 


~ 2 x [K] x [Xe] x exp (~ Ex (X)/T g > | 


+ °"stim, Di^ 


1*3 * (K 


<- V 


eq A, Di 


. x 12) 


0. 25 x [Kr * exp (-V^ ^(M/Tq) 


^stim^Ex 

Di 


= (A e /2c) x 2 x R^ x x 


Di 


Ex ’ 
Di 


( 11 ) 

( 12 ) 


where X and v are the transition wavelength and frequency, c is the 
speed of light, R Ex I s the internuclear separation on the excimer or 
dimer band corresponding to the wavelength, X, A— is the Einstein A 
coefficient assumed equal to the resonance line at coefficient for both 
the excimer and dimer bands, V{X)s are the potential energies of the 
molecular state measured relating to their dissociation products, and 
K and K are the equilibrium constants for the excimer 

6 C| y A f 6 Cj^ y A. f 

and dimer A states. This notation follows that used in Ref. 21. 
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The computer program Listings in the appendix show how the 
coupled rate equations governing all of the above processes are inte- 
grated numerically using a time step chosen by the user and aLtered, 
if desired, during the course of the computer run. The computer pro- 
gram is written for the PDP10 computer. 

B. Results and Discussion 

1 . Discharge Stability Considerations 

Discharge stability considerations are, of course, of 
paramount important in considering cw operation of a discharge 
pumped excimer laser. There is a wide range of discharge instabili- 
ties that can occur m discharge pumped lasers. Many of these instabil- 
ities, such as the various thermal instabilities, have relatively slow 
growth rates and can be stabilized by simply controlling the gas tempera- 
ture or gas flow as is necessary m any case discussed above. How- 
ever, there is one fairly fast growing instability that is a particular 
problem with discharge pumped excimer lasers and deserves special 
discussion. This is the two-step ionization instability which results 
from the high metastable or resonance level population required m 
these systems. Collisional ionization of these levels dominates the 
ionization of an excimer system pumped over laser threshold and Leads 
to a quadratic increase of the ionization rate with electron density. 

The ionization is thus unstable if the electron loss rate varies only 
linearly with electron density as is the case for attachment or diffusion- 
controlled discharges. Recombination controlled discharges are only 
marginally stable against two-step ionization instabilities. 

There are only four basic ways to achieve excimer laser 
threshold without encountering the above instability: 

1. Establish the threshold input power in the discharge 

m a time that is short compared with the growth 
time of the instability (TEA and Bloom line fast 
discharge devices). 
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2 . 


Operate the discharge at below avalanche E/P 
values where the instability does not occur. 

This requires the use of an external source of 
ionization. 

3. Stabilize the instability with a current limit- 
ing element in the discharge circuit, such as 
a ballast resistor. 

4. Operate the discharge in a regime which de- 
pletes the ionizing species thereby saturating 
the two-step ionization rate. 

Approach No. 1 is a demonstrated successful means of obtain- 
ing discharge pumped laser action in the XeF and KrF excimer systems. 
It is clearly inapplicable to cw operation. Approach No. 2, using an 
electron beam as the external source of ionization, is an approach cur- 
rently being pursued elsewhere for high average power operation on the 
rare gas monohalide excimer lasers. It is considered marginal at best 
for cw operation because of thermal loading of the e -beam foil window. 

Approach No. 3, of course, has been used for decades to 
stabilize gas discharges. It is the approach chosen for this program 
for stabilizing the XeF capillary discharge and is the anticipated 
approach {through the use of distributive resistive electrode) to stabiliz- 
ing the transverse XeF discharges in the forthcoming phase of the 
program. 

Approach No. 4 can be used for the alkali excimer system 
because required alkali concentration (which is the ionized species) is 
enough to permit ionization saturation effects to occur. 

Approaches 3 and 4 are the chosen mechanism for stabilizing 
the two-step ionization instability in this program. The numerical 
modeling results, in fact, dramatically reveal the stabilizing effects 
of both of these mechanisms. In the case of ballast resistor stabiliza- 
tion a simple model is presented below which exposes more clearly 
the basic physical requirements for this method of stabilizing the two- 
step ionization instability. 
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The two-step ionization instability in an excimer laser is 
described by the following coupled rate equations for electrons and the 
metastable (or resonance) level population: 


d 

dt 

| ^ 

n = R. n*' n - Bn 
e i e e 

(13) 

_d_ 

dt 

n ~ = R m n e “ n / T ’ 

(14) 


where n^ and n are the metastable (or resonance) level concentration 
e 

and electron concentration, respectively, R. is the rate constant (per 
n") for ionization of the metastable, R^ is the production rate constant 
for metastables, (3 is the linear loss rate of electrons, and is the 
effective lifetime of the me ta stable. 

The steady- state solution to these equations is 


n = S/R R t 
e i m 

o 


(15) 


and 


n = R t n 
o m e 


(16) 


This solution is unstable as shown by linearizing eqs. (13) and (14) 
with respect to the perturbed solution: 


n _ = n + n (t) exp (at) 

C C 6 

o 


(17) 


and 


n 


'I' -v 

n + n 
o 


(t) exp (at) . 


(18) 


This procedure gives 


a = l/2 x l/r x 


-1 ± yi + 4|3 t 


(19) 


Thus the steady state solutions (15) and (16) are unstable with a growth 
rate given by the positive root for a . However, if there is a ballast 
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resistor in the discharge circuit, the rate constants R and R. wiLI be 

m i 

decreasing functions of n^, and one must include this dependence in 
the linearized stability analysis. The most sensitive dependence will 
be on R since it contains the most strongly varying Boltzmann factor 

1XX ] 

(exp(-E m /Te)) compared with expf-fE^ - E m )/ Te ) for R^ where E^ and 
E^ are ionization and metastable level energies, respectively. Thus, 
in the perturbation analysis of eqs. (13) and (14) we also write 

d R 


R 

m 


= R 

m 


+ 


m 

d n 


n 

e 


( 20 ) 


o e 


After collecting terms linear in n as before, the solution for will 

e 

now have an added term: 


a = 1/2 

To determine dR /d n we now make use of the fact that under 
m e 

conditions required for pumping excimer lasers, the dominant electron 
energy loss is due to the inelastic collisional pumping of the meta- 
stable level. Thus the energy conservation relation for electrons can 
be written approximately as 


1/t -1 ± \ 1 + 4pr 


(dR /dn ) 
m e 

Rm 


n 


( 21 ) 


J x E = n R E , (22) 

e m m ' ’ 

where J and E are the current density and field in the discharge. If 
there is a ballast resistor m the circuit this equation can be written 

eN V ,(V - IR)/f = n R E (23) 

e d ' e m m ' ' 

where V is the maintained power supply voltage, I is the current, R is 
the ballast resistance, f is the discharge gap, and is the drift 
velocity. 

Using I = n V , eA where A is the discharge area and e the elec- 

0 Cl 

tron charge we can write 
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R 

rr 
n 

e 

(24) 

Using eq. (21) we see that stability is achieved if Rl/(V - IR) >1, i. e. , 
the drop across the ballast resistor must be equal to or greater than 
the discharge voltage. The onset of discharge stability when this con- 
dition is met is revealed in the numerical modeling for the XeF system 
as well as m the observed current voltage waveforms m the capillary 
discharge pumped XeF experiments. 

2. Modeling Results for the Capillary Discharge Pumped 

XeF System 

A typical result of the numerical integration of the 
modeling equations governing discharge pumping of the XeF excimer 
system ts shown in Figure III— 1 . The curves are discharge current and 
voltage and the excimer level population density versus time. In the 
inset of the figure a typical experimental observation of the same 
temporal profile is shown for comparison. 

As can be seen, the qualitative features of experimentally 
observed waveforms are well reproduced by the theory. The time 
scale over which these features occur is shorter for the theoretical 
case because of the high currents reached in modeling. In both the 
experimental and theoretical waveforms three distinct phases of the 
temporal development are revealed. An initial transient phase, an 
intermediate steady state phase and a final terminated-fluore scence 
phase . 

The initial phase in which the excimer fluorescence reaches its 
peak value is the transient period during which the discharge field and 
electron temperature attain their peak values. These values are 
limited primarily by the discharge circuit inductance. This is the dis- 
charge phase during which XeF excimer laser action is being achieved 
in the TEA and bloomline fast discharge devices. While it may be 
possible to achieve laser action in the capillary discharge during this 


d R 


m 


d n 


V 
c_ 

i E 


m 


RA V d e 


R RAV,e 
m d 

o 

V - IR 
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RELATIVE UNITS 


5117-21 



t, /xsec 


Figure III- 1 . Capillary discharge pimped XeF excimer results: theory versus experiment. 




phase through optimization of the discharge risetime, what happens 
during this transient phase is immaterial for cw operation, and thus 
we give this possibility no further consideration. 

The intermediate steady-state phase of the discharge is the 
desired cw operating regime. Although not shown in the plot, the 
modeling results confirm that the gas temperature rise has been 
stabilized by thermal conduction to the wall during this phase. Also, 
as can be seen, the discharge is stable during this period, the two- 
step ionization instability having been stabilized by the 25 kf2 ballast 
resistor. The magnitude of the excimer fluorescence and discharge 
current during this and the preceding transient phase of the discharge 
are lower in the experiment than those predicted by theory by factors 
of ten. The fluorescence during this phase would have to be increased 
by about a factor of ten to give a practical laser gam coefficient of 
~1 % per centimeter. 

For cw operation a more serious problem than the low value of 
observed fluorescence is the observed and predicted termination of 
fluorescence. This occurs in both the experimental and theoretical 
cases long before the discharge capacitor is drained and has b een 
traced to the consumption of fluorine by the XeF excimer level pump- 
ing reaction itself (reaction G18X). We have identified, therefore, a 
fundamental constraint to cw operation of an XeF excimer laser. The 
basic requirements for overcoming the fluorine consumption problem 
can be appreciated by examining the rate equation for F ? . Replenish- 
ment of at the required rate by gas flow is not practical m a capil- 
lary tube so the rate equation is written for a static gas. 


d[iy| 

dt 


and 


r d!ss r X G ^ - r re= M M = ° 


21 -3 - 1 

4 x 10 cm sec (for 0.2%/cm gain) 


r 1. 5 x 10 ^ cm^ sec 

rec 


(25) 
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To achieve a minimum practical laser gain of 0. 2%/cm the 

21 -3 -1 

pumping rate (consumption rate) must be >4 x 10 cm sec . For 
molecular recombination to stabilize the F 9 concentration at~[F 9 ] ~ 

10 cm where stable discharge operation is presently being achieved 
(higher F^ concentration will lead to unpractically high discharge volt- 
ages to overcome the higher attachment loss), a 50% dissociation frac- 
tion and a third body concentration of ~3 x 10 cm is required. In. 
principle, the requirement of a high buffer gas pressure is not incom- 
patible with the operation of cw capillary discharge pumping of this 
excimer system since neither thermal conduction nor the two-step 
ionization instability is dependent on buffer gas pressure. Experi- 
mentally, however, attempts to elevate the buffer gas (helium) pres- 
sure above -100 Torr have led to filamentmg of the discharge. 

The F^ consumption problem and consequent requirement of 
operating at 1 to 10 atm of buffer gas pressure at present is the most 
outstanding problem for cw capillary discharge pumping of the XeF 
excimer system. Conclusive evaluation of this problem is expected to 
be completed in the early phase of the second year's contract period 
for this program. 

3. Modeling Results for the Capillary Discharge Pumped 

KXe/K ^ Excimer System 

The results of three computer runs of the capillary dis- 
charge pumped KXe/K^ model are presented in Figures III-2 and 
III-3. The series inductance, L, is assumed equal to zero in these 
runs. 

The two results presented in Figure III-l are for a discharge 

through a 10 atm K-Xe mixture at two different temperatures. As is 

the case for optical pumping, the lower temperature result shows most 

o 

of the gain occurring on the excimer band (8500 A), while the higher 
temperature case shows the highest gain occurring on the dimer band. 

In both cases the model predicts that a fairly high E/N value of ~10 ^ 
V/(cm x [K] cm ^) is required to give practical gain coefficients of 
~1% cm. 
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EXe] = 1 X 10 20 cm -3 

C =00 


4u v. = 0 



a. b 

Figure 111 - 2 ,. Theoretical modeling results for avalanche sustained 
discharge pumping of the K-Xe/K 2 system, (a) Low 
temperature results, (b) High temperature results. 

The results m Figure HI- 2 are for an infinite storage capaci- 
tance and no ballast resistor. These results illustrate an important 
feature of the discharge kinetics. Unlike the rare gas monohalides the 
alkali rare gas discharge in the modeling appears to be intrinsically 
stable against ionization instabilities. This stability is apparently due 
to two causes: first, the ground state potassium population is depleted 
during the discharge, resulting in a saturation of the ionization and, to 
a lesser extent, the discharge is recombination stabilized rather than 
attachment stabilized. (Ambipolar diffusion is negligible at these 
pressures. 
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t,^i.sec 

Figure III — 3 . Capillary pumping: low pressure case. 


The low pressure modeling results in Figure III-2 show gain 
only on the dimer band, of course, and the magnitude of the gam coeffi- 
cient is an order of magnitude lower than that of the high pressure case 
under otherwise similar conditions because of the reduced contribution 
of the three-body resonance level associations reaction (reaction GPD). 
Also, in the modeling, a ballast resistor is necessary to stabilize the 
ionization instability, probably because the low pressure discharge is 
wall stabilized (linear loss rate of electrons) rather than recombina- 
tion stabilized as in high pressure case. However, it is not yet cer- 
tain that a divergence of the numerical integration procedure is not 
involved in this case. 

In all three cases modeled gas temperature is adequately con- 
trolled by thermal conduction to the capillary tube wall, as originally 
anticipated m the preliminary evaluation phase of the program. 
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IV. 


LABORATORY APPARATUS 


The objective of this task is to design and fabricate experimental 
apparatus for the purpose of evaluating the potential laser systems 
recommended from our theoretical studies. As discussed in Sections 
II and III, this analytical work identified the potassium-xenon (KXe), 
or argon, excimer, potassium dimer (K.,), and xenon fluoride (XeF) as 
promising laser species to meet the NASA high-power laser require- 
ments for continuous operation in the visible wavelength spectral 
region. Electrical discharge excitation in a small bore capillary tube 
was also predicted to be the best pumping technique to use for pre- 
liminary experimental evaluation of these systems. 

Both the alkali excimer-dimer and the rare gas halides pose 
significant materials problems which must be considered in the design 
of apparatus. Potassium is chemically reactive with many materials 
and must be heated to ~300°C to achieve sufficient density for laser 
oscillation. Molecular fluorine is difficult to handle because it is one 
of the most powerful oxidizing agents known. With these factors in 
mind the gas handling apparatus shown in Figure IV- 1 was designed and 
fabricated for use with the XeF laser system. This apparatus is con- 
structed primarily of stainless steel tubing and incorporates stainless 
steel packless valves and fluorine compatible pressure gauges. The 
entire system was passivated with fluorine-helium mixtures and pure 
fluorine prior to use for discharge measurements. Nitrogen trifluoride 
(NF^) presents less of a materials problem and was also used in this 
apparatus as a fluorine source for reaction with excited xenon atoms. 

The discharge tube and electrical pulse circuit are shown photo- 
graphically in Figure IV- 1 and schematically in Figure IV-2. Stain- 
less steel cylindrical cavity electrodes as shown in these figures were 
found to be more satisfactory than simple tunsten pin electrodes which 
were used in our first XeF discharge experiments. The discharge 
tube was situated within an optical cavity as shown in Figure IV- 1. 
Alignment of the optical cavity was achieved with the use of a He-Ne 
laser and conventional alignment techniques. 
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Figure IV- 1. Rare-gas halide capillary discharge laser 
apparatus. 

TO GAS HANDLING 
APPARATUS 



Figure IV-2. Schematic diagram of XeF capillary discharge tube and 
pulse circuit. 
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We have also designed and constructed an apparatus for 
conducting discharge experiments with potassium-rare gas mixtures at 
temperatures between 250 to 300°C. A schematic drawing of the high 
temperature apparatus is shown in Figure IV-3. Separate ovens are 
used for the discharge tube and potassium reservoir in order to elimi- 
nate condensation effects in the discharge tube. Homgeneous mixtures 

3 

of potassium vapor and rare gases are prepared in a 100 cm stain- 
less steel cylinder. Partial pressures of potassium are controlled by 
the temperature of this reservoir which contains enough potassium to 
produce satuated vapor pressures. The small volume discharge tube 
is then filled by gas expansion from this chamber: the precise potas- 
sium dimer (K ) absorption at 6328 A and resonance line absorption 
^ o o 

by atomic potassium at 7665 A and 7699 A. A photograph of this 
apparatus is shown in Figure IV-4. The remainder of the electrical 
and gas handling apparatus for the potassium experiments is shown in 
Figure IV- 1. 

All of the discharge experiments conducted on this program 
have been done with the apparatus shown in Figures IV - 1 and IV-2 and 
will be discussed in the following section. Absolute fluorescence mea- 
surements were conducted with a variety of gas mixtures and electri- 
cal parameters in order to evaluate discharge characteristics and 
laser potentialities. The optical configuration used for this purpose 
is shown schematically in Figure IV-5. A calibrated photodetector and 
optical filter were used to measure the side light emission from a 
well-defined volume in the discharge tube. These measurements can 
then be used to determine excited state populations of XeF and the 
maximum optical gain coefficient, a, with the assumption that the 
lower laser level is unoccupied. The relationships used to obtain this 
information are given below. 
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Figure IV-3. Schematic diagram of potassium rare gas discharge apparatus. 
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Figure IV-4. Photograph of potassium rare gas apparatus in final stages of 
fabrication. 
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CAPILLARY TUBE 


un 

o 


7.6 cm 



SHIELD 


0.3 cm SLIT 



CORNING CS-7-39 FILTER 

UDT CALIBRATED PHOTODIODE 
(PIN-10 UV; 1.14 cm diam) 

TO BIAS BOX AND 
OSCILLOSCOPE 


DETECTOR ENCLOSURE 


GEOMETRICAL COLLECTION FRACTION e F = 1.37 x 10~ 3 


Figure IV-5. Experimental arrangement for absolute fluorescence measurements. 
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where 


P^, = total optical power emitted from volume, V 
Ip = photodetector current 

F = geomtrical collection fraction (1.37 x 10 ) 

Sp = detector sensitivity (10 W/A at 3450 R.) 

T = average transmission of filter (0.464 from 
340 A to 3500 A) 

A = Einstein spontaneous emis sion coefficient 
(2.0 x 10? sec"^) 

AX = full width at half maximum of XeF^ emission 
band (105 A). 

Use of the appropriate constants with the. expression above shows that 

13 - 3 

an upper state XeF concentration of 3.76 x 10 cm is required to 

- 3 - 1 

give a small signal peak gam coefficient of 1 x 10 cm 

Spectral emission profiles were measured with a l/4 m mono- 
chromator equipped with an optical multichannel analyzer detection sys- 
tem. The fluorescence spectrum of XeF obtained with this apparatus is 
shown in Figure IV-6. An NF^XerAr = 1:1:10 mixture at a total pres- 
sure of 30 Torr was used in the quartz capillary discharge tube to 
obtain these data. 

4950-4 


Figure IV-6. 

Spectral profile of 
XeF excimer fluores- 
cence from capillary 
discharge tube. 

333 343 353 363 373 

X r nm 
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V. 


EXPERIMENTAL RESULTS 


A Introduction 

This section contains the results of Task IY — Molecular 
System Evaluation. The objective of this task is to experimentally 
evaluate the two excimer systems recommended m Task II (K^/KXe 
and XeF) with the apparatus designed under Task III (see Section IV). 
We describe measurements taken on the XeF capillary discharge 
apparatus. Operation of the K^/KXe metal vapor capillary was delayed 
because of leak problems and no data have yet been obtained. It is 
anticipated that the metal vapor experiment will be operational early 
in the next contract period. 

The primary measurements carried out were tests for gain 
and laser action, the determination of operating conditions, limiting 
processes, and scaling information, These latter two points are dis- 
cussed in Section III. Since laser action was not achieved on XeF, 
additional tests were made to determine why it did not perform as pre- 
dicted, with the objective of gaming new knowledge of the XeF system 
in order to improve the analysis of Task II (Section III). 

Data obtained from the XeF apparatus included discharge volt- 
age, discharge current, and absolute fluorescence intensity measured 
as a function of gas mixture, pressure, and external circuit param- 
eters (e.g., capacitor and ballast resistor size). Upper state XeF 
populations and small signal gain were calculated from the fluorescence 
data according to the procedures outlined m Section IV. 

Experiments were performed with two different size capillaries: 
(1) a 15 cm long x 0. 3 cm diameter tube, and (2) a 1. 6 cm long by 
0. 1 cm diameter tube. Tube 1 was used in the initial experiments and 
was designed to provide a gam path long enough for laser action and a 
bore size large enough to permit easy mirror alignment. Tube 2 was 
used m the second phase of the experiment primarily to obtain a better 
understanding of the system kinetics. This tube was designed to pro- 
vide more efficient gas cooling and higher e/P than tube 1. 
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Several unsuccessful attempts at pulsed laser action were made 

with Tube 1 under conditions where the fluorescence measurements 

-3 - 1 

indicated peak gams of 8 x 10 cm . The modeling studies indicated 
that the gas temperature was too high and E/P too low to produce an 
inversion. The high fluorescence yields observed were apparently 
dominated by thermal excitation. 

Tube 2 was then used to obtain additional data under more favor- 
able operating conditions (lower temperature, higher E/P) which were 
then used as input to the modeling study. The pulsed and cw perform- 
ance of Tube 2 was not significantly better than Tube 1. Pulsed gams 

- 3 - 1 

were typically 0. 8 x 10 cm for~2 psec durations. Continuous wave 

-5 - 1 

gains were typically 1x10 cm and limited in duration only by the 
thyratron turnoff characteristics and fluorine consumption as discussed 
in Section III. 

B . Experiments With Tube 1 

The initial experiments with Tube 1, the 15 cm long x 0. 3 cm 
i. d. tube, were carried out with a 50 to 200 kf2 resistor in series with 
the discharge capacitor to investigate the quasi-cw behavior of the sys- 
tem. Eater experiments were performed with R = 0 to examine the 
short pulse behavior. A summary of the experimental conditions and 
main results of these studies is shown in Table V-l. 

1 . Quasi-CW Experiments With Tube 1 

Both He and Ar were tried as diluents and NF^ was 
employed as the fluorine source. The He mixtures typically exhibited 
slightly better ('"•'20%) fluorescence intensities than comparable Ar mix- 
tures. A typical mixture ratio used m these experiments was 1:1:2/ 
NF^-Xe-He. Uniform quasi-cw discharges of 100 psec duraction could 

be readily obtained at total pressures 10 Torr, at current densities of 
2 

1 to 10 A/cm and E/P values of ~100 V/cm Torr, At pressures 
higher than ~10 Torr, the discharge would detach from the capillary 
wall and run a twisting course through the tube. The fluorescence 
intensity would be decreased by a factor of '"'2 when in this mode. 
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Table V-l. Operating Conditions and Main Results of 
15 cm Long x 0.3 cm i.d. Tube Experiments 

Mixtures NFg.Xe-He = 0.1 - 1 5: 1 - 3 1.67-10 

NF 3 ‘Xe:Ar= 1 - 1.5: 1 -3:10 

NF 3 .Xe =1:4-9 

Pressu res 0. 5 — 1 0 T orr 

Apparatus C = 260 and 780 pF 
R = 0 to 200 
20 0 kV charging voltage 

Long Pulse Experiments 

XeF* of order 10^ — 10^cm — ^ (a~ 1 0 ' 5 cm“^) 
Fluorescence Pulses ~ 50 + /usee long 
Current Densities J ~ 1 — '1 0 A cm^ 

Steady State E/P ~ 100 V/cm — Torr 

Short Pulse Experiments (R — 0) 

XeF* Peak Values of 1 to 3 x 10^ cm - ' ^ 

Apparent Gains of a - 2.6 — 8 0 x 1 0 — ^ cm - ^ 

Fluorescence Pulses ~ 0.5 psec long 

Current Pulses ~ 0 2 /usee long, J = 100 to 4000 A cm^ 

E/P Low (Unmeasured) 

Laser Action Tried With 1 to 4/NF 3 — Xe Mixture 
Attempt Unsuccessful Because of Gas Heating 

T1938 
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2 . 


Short Pulse Experiments with Tube 1 


a, Fluorescence Optimization — As the ballast 
resistor (R) is decreased the current pulse shortens, the initial cur- 
rent rises and the fluorescence changes from a quasi-cw to a pulsed 
character. The peak fluorescence increases with increasing current 
all the way to the point where R = 0 and the current is limited primarily 
by the circuit inductance. The increase m fluorescence obtained by 
going from R = 100 or 200 to R = 0 was typically a factor of ~200. 

Further increases m the fluorescence intensity were obtained 
by varying the mixture ratio and gas pressure. It was found that the 
fluorescence increased with increasing NF, and decreasing He diluent. 

t* 

Ultimately, conditions were found where peak XeF populations of 
14 -3 -3 -1 

3x10 cm and gains of 8 x 10 cm were observed. This corre- 
sponds to a single pass gain of 12% in our 15 cm tube. 

Figure V-2 shows fluorescence and current traces for this 
maximum gain condition in the 3 mm tube. The fluorescence from 
this l^/NF^-Xe mixture at 7 Torr has an FWHM of ~300 nsec. The 

current pulse is also '■"300 nsec wide and reaches a peak value of 
/ 2 

'"'700 A/cm . The discharge voltage was unmeasurable because of 
noise, but the E/P was probably 50 V/cm-Torr. In Figure V-2 the 
fluorescence trace is displaced 2 divisions to the right relative to the 
current trace. 

b. Laser Action Experiments — Laser action was 

attempted with this l^/NF^-Xe mixture at pressures of 1 to 10 Torr. 

The optical cavity consisted of two 50 cm radius mirrors spaced 40 cm 

apart. The mirrors were made of quartz and coated with multilayer 

dielectric films to give high reflectivity and 0. 1% transmission near 
o 

3500 A. This stable cavity was carefully aligned with an external He- 
Ne beam. A photomultiplier equipped with a narrow band 3500 X filter 
was placed '"’40 cm from the output window in the output beam path. 

The apparatus was carefully apertured so that the phototube saw only 
the light emitted through the output mirror. The detector system was 
sensitive enough to easily see the nonlasing XeF'-! 5 spontaneous emis- 
sion obtained when the rear laser mirror was blocked. 
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1-4/NF 3 -Xe 7.0 Torr 

780 pF AT 20 kV 
NO DISCHARGE RESISTOR 
TIME 0.5/isec/div 

FLUORESCENCE (TOP) 1.03 x 10 14 XeF* cm _3 /div 

XeF* pEAK = 3.0x 10 14 cm -3 

3 i 

a p ~ 8 x 10 cm 
CURRENT (BOTTOM) 100 A/div 

'peak- 700 A/cm2 

Figure V-2. Optimi ed fluorescence for 
3 mm tube experiments. 


The onset of laser action with the above system would have been 
observable as a drastic (several orders of magnitude) change m the 
photomultiplier signal level. By comparing the detector signals 
obtained with the rear laser mirror blocked and unblocked, even 
extremely weak laseing action should have been observable. 

No change in the phototube signal was observed despite pres- 
sure and current variations as well as mirror excursions performed 
about the He-Ne alignment condition. Laser action with a 1. 5:3.0: 
10/NF^-Xe-He mixtures was also tried with a similar lack of success. 
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As was explained in Section III, the lack of laser action is 
probably due to a high gas temperature and low E/P brought on by the 
high current densities required to give apparent gains of practical mag- 
nitude. The high fluorescence yield is not due to selective excitation of 
XeF*, but instead is the result of thermal pumping which also populates 
the XeF ground state. ~ 

c . Experiments With the 1, 6 cm. Long x 0. 1 cm i. d. 

Tube (Tube 2) — Quasi-cw as well as short pulse 

operation was also examined m the 1 mm tube. Current densities were 

2 

usually kept below 50 to 7 5 A/cm to prevent significant gas heating. 

A summary of the experimental conditions and main results of the 1 mm 
tube studies are shown m Table V-2. 

1 . Quasi-CW Experiments With Tube 2 

a. Measurements on NF ^ and F , Mixtures — The 

first experiments with the smaller tube were directed toward quasi-cw 
operation. Both F, and NF were used as fluorine sources and He was 

^ J 

used as a diluent. A typical mixture inve stigated was 0. 2:1: l/NF^fFM- 
Xe-He. Data for this mixture ratio are shown in Figure V-3 {with NF^) 
and Figure V-4 (with F ). Total anode-cathode voltage drop, XeF* 
fluorescence intensity and total discharge current are shown as func- 
tions of total pressure. Operating pressures were usually limited to 
below 30 to 40 Torr by the discharge not breaking down. 

The NF, data (Figure V-3) show that after an initial 50 to 100 

^ 12 -3 

psec high fluorescence period with XeF* = 1. 5 x 10 cm , the dis- 
charge voltage and current reach a steady state and produce a quasi-cw 
XeF* population for ~0. 5 msec. The termination of the fluorescence 
is caused by the turnoff of the thyratron and is not due to kinetics 
effects. Rough calculations indicate that the fluorescence could remain 
at these levels for several milliseconds before NF^ consumption effects 
become significant. 
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Table V -2. Operating Conditions and Main Results of 
1. 6 cm long x 0. 1 cm i. d. Tube Experiments 


Mixtures 

NF 3 :Xe-He = 0 1 to 1:1:1 - 2 


F 2 :Xe:He = 0.2*1 :1 

Pressures 

0.5 - 40 Torr 

Apparatus 

C = 260 and 780 pF 
R = 5 to 500 Kf2 
20 kV Charging Voltage 

Long Pulse Experiments 


Steady State XeF* of 2 — 4x10^ cm 2 Observed for 
0 5 msec Durations in NFg and F 2 Mixtures 


Steady State Gains = 05x1 0 — ^ — 10x1 0 — 5 cm~^ 
J = 1 - 5 A cm 2 


E/P = 60 to 250 V/cm - Torr 

Short Pulse Experiments (R > 5 K£2) 


Peak XeF* of ~ 3 x 10^ 2 cm - ' ^ (a ~ 0.8 x 10 — 2 cm - ^ ) 
Fluorescence Pulses 1 to 20 psec Long 
J = 5 — 50 A cm 2 


Discharge Runs at High E/P (400 to 700 V/cm — Torr) 
for ~ 5 psec at Start 


Discharge Then Makes Transition to Lower E/P (80 to 300 V/cm — Torr) 
Maximum Fluorescence Occurs During High E/P Initial Phase 


T1939 
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0.2 1.1/NFg- Xe- He 
FLUORESCENCE: 4.71 x 10 11 cm -3 /div 



Figure V-3(b). Continued. 



0 2.1 :1/F 2 — Xe — He 

FLUORESCENCE = 1 ;1 2 x 10 12 cm -3 /div 



Figure V-4(b). Continued. 
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The steady state fluorescence scales approximately linearly 

11-3 

as the total pressure from a value of 1. 2 x 10 cm at 7. 5 Torr to 
11 -3 

3.7 x 10 cm at 23 Torr. The corresponding quasi-cw gains are 
-5 -1 -5 -1 

0. 3 x 10 cm and 1.0 x 10 cm , respectively. Steady state E/P 
falls from 250 V/cm-Torr at 7. 5 Torr to 135 V/cm-Torr at, 23 Torr, 
Current densities are about 2 A/cm for all three cases shown. 

It should be noted that the current does not decay as fast as is 
indicated m Figure V-3. About half of the decay is caused by the 
0. 06%/jxsec droop (36% in 600 psec) of the Pearson 2100 Rogowski coil 
used in these measurements. This spurious decay is quite evident 
at the end of the discharge where the current trace is observed to 
undershoot the baseline. 

The F^ data (Figure V-4) show similar time and pressure 
variations. After an initial transient, a 0. 5 to 0.7 msec quasi-cw 
fluorescence level is observed during a period of nearly constant dis- 
charge current and voltage. The steady state fluorescence again 

11 -3 

scales approximately linearly with pressure from 1x10 cm at 

15 Torr (0..27 x 10 ^ cm * gain) up to 4. 5 x 10'*' 1 cm ^ at 40 Torr 
-5 - 1 . 

(1.2 x 10 cm gain). The mixtures containing F^ thus produce 
about one-half as much fluorescence as corresponding ones containing 
NF^. This may be a result of the visually poorer discharge uniformity 
obtained with the F^ mixtures. Whereas the NF^ mixtures were uni- 
form across the entire tube bore, the mixtures were noticeably 
brighter at the tube center than they were near the tube wall. The 
steady state E/P was also lower in the F^ mixtures, decreasing from 
80 V/cm-Torr at 15 Torr to 60 V/cm-Torr at 40 Torr. 

b. Fluorine Atom Recombination Experiment — As 

discussed m Section III, a significant factor limiting cw operation is 
consumption of the fluorine source by the reaction Xe* + FHNF-) -**- 
XeF + F(NF/. One way to circumvent this problem is to add a large 
amount of diluent to the mixture to enhance the three body recombina- 
tion reaction F + F + He — + He. To see if this recombination reac- 
tion would have any effect upon the fluorescence level observed during 
quasi-cw operation, a 0 . 2: 1 :40 /F^-Xe-He mixture was studied at 
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pressures up to 200 Torr. A 0. 2:1:40 mixture at 200 Torr has the same 
and Xe concentrations as a 10.7 Torr 0. 2:1:1 mixture, but 40 times 
more He to act as a third body. 

Figure V-5 shows the 200 Torr operation of this mixture. 

The gas broke down easily, but the discharge was filamentary in nature 
with only about one-half of the bore filled with a bright discharge 
directly down the tube center. The nonumformity in the case was much 
worse than m the previous 0.2: 1:1 mixture. As before, the discharge 
lasted for ''■'0.7 msec at constant voltage and approximately constant 
current until the thyratron shutoff. e/P was 6.25 V/cm-Torr (258 V/ 
cm-Torr Xe), and assuming that the current and XeF-^ population were 

2 

uniformly distributed throughout the tube, these had values of 2. 5 A/cm 
11-3 

and 2. 35 x 10 cm , respectively. 

This XeF* population is about a factor of two greater than that 
observed m 0. 2:1:1 mixtures at 10 to 15 Torr, indicating that the added 
He may be having a positive effect. However, the filamentary nature 
of the discharge makes these results suspect. The '‘filament" could 
be an "arc" with a high current density and a thermally excited popula- 
tion, similar to the high current experiments with the 3 mm tube. 
Further studies are required before a definite conclusion can be made 
on the efficacy of the recombination approach. 

2. Short Pulse Experiments With Tube 2 

a. Current Density Limitations — As the current 

is increased in the 1 mm tube (Tube 2), the initial fluorescence intensi- 
ty also increases. The current cannot be increased indefinitely, how- 
ever, since eventually gas heating becomes excessive and any hope of 
obtaining a XeF inversion is destroyed by a lowered E/P and thermal 
filling of the XeF ground state. We therefore conducted several 
experiments in an effort to determine an upper limit on the current 
density below which high E/P and good XeF* fluorescence could be 

2 

obtained simultaneously. This limit is found to be ~50 to 75 A/cm 
for the mixtures studied. 
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Figure V-5. Operation at 200 Torr with high helium concen- 
tration. 


67 


At current densities below this limit, high E/P is observed for 
~5 (isec at the beginning of the current pulse. After this initial high 
E/P phase, the discharge then makes a transition to lower E/P and 
continues to operate there until the thyratron shuts off. The high E/P 
phase is characterized by E/Ps in the range 400 to 700 V/cm-Torr, 
while 80 to 300 V/cm-Torr is typically observed during the low E/P 
phase. The maximum XeF* fluorescence is observed during the low 
E/P phase. The maximum XeF* fluorescence is observed during the 
high e/P phase. The exact values of E/P obtained depend upon mix- 
ture ratio and pressure. E/P generally increases with decreasing 
pressure and increasing NF^IF^). 

This dual mode of operation is illustrated by photographs 1, 2, 
and 3 of Figure V-6. This figure shows discharge voltage and current 
for decreasing values of ballast resistor. An initial ~5 psec voltage 
hump (the high E/P phase) is followed by a relatively flat voltage trace 

(the low e/P phase). Despite the fact that the current density varies 
2 2 

from ~10 A/cm in photograph 1 to ~40 A/cm in photograph 3, all 

three cases exhibit e/Ps of ~400 V/cm-Torr during the first phase 

and ‘~~T25 V/cm-Torr during the second. 

This collapse of the initial high E/P is primarily a result of 

consumption of NF^F^) as discussed in Section III. Since this effect 

is proportional to current density, one would expect the high E/P phase 

to shorten drastically at high current densities. This, in fact, is 

2 

observed. At current densities above 50 to 75 A/cm most mixtures 
exhibit a high E/P phase of only several hundred nanoseconds duration. 
This effect is illustrated by photographs 4 and 5 of Figure V-6, where 
initial current densities are ~150 and ~350 A/cm , respectively. The 
E/P m these mixtures falls almost immediately to low values of 
~80 V/cm-Torr. At current density of this level, excessive gas heat- 
ing will also have a detrimental effect upon the XeF* inversion. 
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0 1 A/div 


R = TOO kS2 


5 kV/div 

0.2:1:1/NF 3 -Xe-He 
1 5 Torr 
1.0 mm BORE 

20 kV CHARGING VOLTAGE 

5 kV/div 


R = 50 kft 


0.2 A/div 


5 kV/div 


R = 25 kfl 


0 2 A/div 



u 


5 /usec/div 


Figure V-6(a). V-I data for NF^ mixtures illustrating cur- 
rent density limit. 
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5 kV/div 


R = 1 0 k£2 


1 A/div 


, - Xe - He 


RGING VOLTAGE 


5 kV/div 


R = 5 kS2 


1 A/div 


Figure 



b. Discharge Stability — In Section III a stability- 
condition was derived which indicated that the discharge will be stable 
as long as the voltage drop across the ballast resistor is greater than 
the voltage drop across the discharge. The data shown m Figure V-7 
tend to support this stability condition. 

This series of oscilloscope traces shows the V-I characteristic 
of a 0. 2: 1: l/NF^-Xe-He mixture as a function of pressure for a fixed 
charging voltage (20 kV) and fixed ballast resistor (50 kl2). In photo- 
graphs 1 through 3, the peak discharge voltage is well over half of the 
total 20 kV on the discharge capacitor. It is 15 kV, 13 kV, and 11 kV 
in photographs 1, 2, and 3, respectively. The ballast resistor thus 
has ~5 kV, ~7 kV, and ~9 kV across it m the corresponding cases. 

The discharge should be unstable, and indeed the discharge current 
does exhibit a factor of 2 type fluctuations in the first 10 psec of each 
trace . 

In photographs 4 and 5, however, these fluctuations are absent. 
In these cases the peak discharge voltages are 8 kV and 5 kV corre- 
sponding to ballast resistor voltage drops of ~12 kV and ~15 kV. These 
voltages satisfy the stability criterion and indeed the V-I characteris- 
tics of 4 and 5 do come to a steady state more rapidly with less severe 
fluctuations than do cases 1, 2, and 3. 

c. Fluorescence Measurements — With the 50 to 

2 

75 A/cm limit in mind, variations m system parameters were carried 
out m order to maximize the XeF* fluorescence. The XeF* population 
was found to be a weak function of current density, pressure, and mix- 
ture ratio. Typically, the peak fluorescence would vary by only fac- 
tors of 2 to 4 for fairly large excursions in operating conditions. The 
use of F^ in place of NF^ also produced only small differences in 
performance. 

Figure V-8 illustrates how the XeF* fluorescence varies with 

pressure for a typical NF -Xe-He mixture. The optimum pressure 

a 13 _3 

for this mixture is ~10 Torr where XeF* = 3.3 x 10 cm (l. e. , gam 
-3-1 

= 0.9 x 10 cm ) is obtained. The fluorescence, which is signifi- 
cant only during the initial part of the current pulse (high E/P phase). 


71 









5117-36 



Figure V-7(b). Continued. 
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20 Torr ALLTRACES- 0 5 1:1/NF 3 - Xe - He 15Torr 
2 jusec/div 


FLUORESCENCE (TOP = 3.3 x 10 13 cm -3 /div 
CURRENT (BOTTOM) = 0 2A/div 



Figure V-8. Pressure dependence of fluorescence from a typical NF^ mixture. 




varies by only about a factor of two over the range 5 to 20 Torr. The 

2 n 

current density varies from ~25 A/cna at 20 Torr to ~65 A/cnr at 
5 Torr. 

Figure V-9 shows oscilloscope traces of the maximum fluores- 
cence obtained for four different X:l:l/NF^-Xe-He mixtures with X = 
0.2, 0.35, 0.5, and 1.0. These mixtures produced the highest fluores- 
cence of any evaluated during this part of the study. In this series, the 
operating pressure has been varied with constant charging voltage 
(20 kV) and a fixed ballast resistor (25 kD) until maximum fluorescence 
was obtained for each mixture. 

Despite the wide variation in NF^ concentration, the optimum 

fluorescence obtained is very nearly the same for all mixtures: XeF* 

13 -3-1 

^3 x 10 corresponding to ''■'0. 8 x 10 cm gain. As indicated 

previously, this lack of a strong dependence upon operating conditions 

appears to be a general characteristic of the XeF* capillary system 

and is not restricted to this particular set of mixtures. 

Figure V-10 shows typical short pulse operation with F_ at two 

2 2 ^ 

different current levels 15 A/cm (top) and 30 A/cm (bottom). Despite 
these relatively modest current densities, the discharge does not 
exhibit the initial high E/P phase observed in NF n discharges. Maxi- 
mum E/P in this mixture is ~200 V/cm-Torr. The fluorescnece ob- 
served is also somewhat lower than in the NF^ mixtures. As in the 
quasi-cw experiments, the F^ mixture discharges were visibly nonuni- 
form with the center being considerably brighter than the regions near 
the walls. This lack of discharge uniformity is undoubtedly linked to 
the relatively poor performance of these mixtures. 
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15 Torr 14Torr 

0.2,1,1/NF 3 — Xe — He 0.351 1/NFg - Xe - He 



ALL TRACES 2 Msec/d iv 


FLUORESCENCE (TOP) = 3 3 x 10 13 cm~ 3 /dtv 
CURRENT (BOTTOM) = 0.2A/div 



10 Torr 10 Torr 


Figure V-9. 


Optimized fluorescence in 1 mm bore tube for several NF 3 mixtures. 





VI. 


PRELIMINARY CONCLUSIONS AND FUTURE WORK PLAN 


A. Preliminary Conclusions 

At this point in the program we are in a position to make some 
preliminary conclusions regarding the expected performance charac- 
teristics of a capillary discharge pumped cw excimer laser. In both 
the XeF and the KXe/K 2 excimer systems the discharge electric field 
required for achieving a practical gain coefficient of ~1% per cm is 
predicted to be high, ~10 to 40 kV/cm. If the field is to be applied 
axially along the capillary tube it will be necessary to use multistage, 
individually ballasted discharge paths as shown m Figure VI- 1 m 
order to limit the required power supply voltage to ''■’10 to 40 kV. This 
type of discharge configuration is used routinely for CC> 2 waveguide 
lasers. The ballast resistors will need to have a voltage drop compar- 
able to or greater than the discharge voltage on each segment so that 
all segments will draw equal current and run in a stable mode. 

In both the XeF and KXe/K^ excimer systems the efficiency is 
predicted theoretically to be near the quantum efficiency ('"’30% for XeF 
and ~90% for KXe/K 2 ). In the case of XeF the optimized experimental 
efficiencies appear to be running about a factor of 10 below quantum 

efficiency. At lasing threshold we expect the input current density to 

2 2 2 
be ~10 A/cm for the XeF system and ~1 A/cm for the KXe/K 2 sys- 
tem. For cw operation both systems will likely have to run at a total 
pressure of between 1 to 10 atm. For adequate gas temperature cool- 
ing the capillary diameter will be required to be ~ 1 to 3 mm (a lower 
figure likely applying to the higher threshold XeF system), with at 
least half of the gas being helium. Flowing liquid coolants, of course, 
will also be necessary to control the temperature of the tube wall. 
Assuming that an efficiency of 10% can be reached for both systems, 
the expected cw power output for near threshold operation should be 
~ 10 to 100 W for the capillary discharge pumped KXe/K 2 system and 
~100 to 1000 W for the capillary discharge pumped XeF system. 




15 ?. 
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Figure VI- 1. Projected capillary discharge configuration for cw exci- 
mer laser. 


B. Future Work Plan 

The capillary discharge technique in use for the present phase 
of the program was chosen as the most expedient means to demonstrate 
and gather basic data for cw operation of an excimer laser system. 

This technique is not dimensionally scalable. Thus, the appropriate 
goals for the future 12 month phase of the program in meeting the over- 
all objectives of the NASA high power laser program are 

1. Identify and show theoretical feasibility of a cw 
pumping technique which can be applied to at 
least one excimer system recommended from 
the current program and which can be scaled 
to produce cw output power sufficient to meet 
NASA's future mission requirements 

2. Demonstrate feasibility of the above concept 
in a laboratory experiment. 

According to the discussion in Section II the only option other 

4 5 2 

than thermal conduction for achieving the 10 to 10 W/cm heat 
removal rate required of a cw excimer laser is convection. For the 
alkali excimer s this is a viable approach to more scalable cw opera- 
tion since residence lengths as long as 10 cm at 1 to 10 atm pressure 
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sonic flow can be tolerated. In the case of the XeF system the 
maximum, residence length required could be as short at 1 mm, which 
would not offer any improved scalability features over the capillary 
discharge. However, there is a particular feature in the XeF excimer 
system pumping kinetics that permits consideration of a new non- 
quilibrium convection and radiative cooling mechanism for achieving 
scalable cw operation of this system at low pressure. 

In the case of the XeF excimer lasers, gas flow can be used in 
a unique way to permit power extraction throughout a dimension >1 cm 
in the flow direction. The method is to utilize a unique feature of the 
kinetics of these excimer systems which will permit one to use low 
pressure gas flow to spatially separate regions of the flow where electron- 
heavy particle thermalization has occurred (downstream region) from the 
region where the electron and gas temperature are still disparate (up- 
stream region). The unique feature is the fast two-body excimer forma- 
tion rate , e. g. , 


Xe* + F- 


XeF + F 


(26) 


which at loiv pressures ~1 to 10 Torr can proceed rapidly enough to 

quickly reach lasing threshold and, m particular, proceeds at a rate 

faster than the electron-heavy particle thermalization rate which 

occurs primarily via vibration-to-translation transfer collisions. The 

-9 -3 -1 

rate for the former process is '"'10 cm sec while the V-T trans- 
fer rate is at least two orders of magnitude smaller. 

17 -3 

Assume a 10:lXe:F_ mixture at a density of 10 cm and 
5 ^ 

flowing at a speed of '"'10 cm/sec. Energy is fed into the electrons at 
the head of the plasma jet and establishes an electron temperature of 
~4 eV (typical of pulsed static XeF£ discharge pumped excimer lasers). 

These electrons quickly establish comparable reservoirs of electronic 

ri 14-3 .. r i 

level energy [Xe 5;; j = 10 cm and vibrational level energy (F1J ~ 

14 -3 L 

10 cm . At these densities and at the above flow velocity, reaction 

(26) will establish an excimer level population sufficient for lasing in 
about 0. 1 mm, i. e. , at essentially the upstream boundary of the plas- 
ma jet, while V-T transfer and thermalization will not be complete 
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until several centimeters downstream {see Figure VI-2). In this 
latter region the gas will become hot (Tg ~ 1 eV) and radiative cooling 
can occur thus giving the added benefit of helping to minimize the heat 
exchanger bulk in a recirculating system. 

Finally, the fluorine consumption problem associated with cw 
low-pressure operation of this system is removed as the gas flow is 
adequate to replenish the fluorine at the required rate. 

In principle, the above concept is similar in many respects to 
many chemical lasers but differs of course m the sense that power 
must be put into the medium by means of a discharge. If, as antici- 
pated, laser power extraction can occur throughout plasma jet volumes 
. , 3 

on the order of 1 cm , a gam path of 50 to 100 cm could provide net 
cw laser power levels in the range that would meet NASA's ultimate 
mission requirements. 

Whether one uses equilibrium or nonequilibrium gas flow to 
establish cw operation of an excimer laser with gas residence length 
of 1 cm, one must In either case first demonstrate stable operation of 
a transverse discharge ina static gas for times ~10 |jsec which roughly 
corresponds to the sonic flow residence time of gas in the discharge. 
This will be done in the forthcoming phase of the program for both the 
KXe/K^ and XeF systems prior to any experiments utilizing gas flow. 

In total, the activities planned for the forthcoming year are 


1 . 


2 . 


3. 


Achieve spontaneous emission brightness levels 
adequate for lasing KXe K or XeF m 

L. 4 


a cw capillary discharge 

a 10 jisec static transverse 
discharge 


Parallel activities 
for first six months 


Construct and achieve cw 
laser action in a long gain- 
path capillary discharge 

Construct and achieve 
several millisecond laser 
action in a fast flow trans- 
verse discharge. 


\ 

Parallel activity 
for second' six months 
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e" + Xe * 

e~ + Xe • 

Xe* 

e“ + F 2 

Xe* + F 2 

e + F 2 * 

F 2 (VIB) + Xe 
F 2 + Xe — 


Xe* + e 
Xe + + e 


F“ + F 
XeF* + F 
e + F 2 (VIB) 


FAST (< 1 fj. sec) 


F 2 + Xe (TRANS) 
XeF g + F 


SLOW (> 10 /a sec) 


5007-2 



*CW OUTPUT POWER = 10 MW/liter 


Figure VI-2. Low pressure XeF excimer plasma jet concept. 
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Te, Tg, eV 




At the end of this second year effort on the program we expect to have 
identified both an excimer molecular system and cw pumping technique 
which is dimensionally scalable to cw power levels on the order of one 
megawatt and to have carried out a small scale demonstration of lasing 
such a system xn the laboratory. 


84 



REFERENCES 


1. C. E. Moore, Atomic Energy Levels , National Bureau of Standards. 

2. H.M. Michaels and F. E. Harris, J. Chem. Phys. 19 , 1464 (1963). 

3. A. J. Palmer, J. Appl. Phys. 41 , 438 (1970). 

4. J, Pascale and J. Vandeplanque, CEA Report, March 1974. 

5. C. York and A. Gallagher, Jt. Inst. Lab. Astrophys. Rpt. 114, 
October 1974. 

6. R. Carbone and M. Litvak, J. Appl. Phys. 39 , 2413 (1968). 

7. R. Hill et al. , Appl. Phys. Lett. 

8. C. Cooper, C. Cobb, and E. Tolnas, J. Mol. Spectrosc. 7_, 223 
(1961). 

9. C. Herzberg, Spectra of Diatomic Molecules (B. Van Nostrand, 

1950). " 

10. J. Ewing and C. Brau, Phys. Rev. Al 2 , 129 (1975). 

11. E. George and C. Rhodes, Appl. Phys. Lett. 23 , 139 (1973). 

12. A.V. Phelps, Jt. Inst. Lab. Astrophys. Rpt. 110, September 1972. 

13. A. von Engle, Ionized Gases (Oxford Clarendon Press, 1965). 

14. L. J. Kieffer, Jt. Inst. Lab. Astrophys. Information Center 
Rpt. No. 13, University of Colorado, September 30, 1973. 

15. S. C. Brown, Basic Data of Plasma Physics (The MIT Press, 

1967. 

16. P. Mahadevan and R. Hofland, Aerospace Corporation reprint 
(1976). 

17. C.A. Beam and J. J. Ewing, Semiannual Tech. Rpt., Contract 
N00014-75-C-0063, Avco Everett Research Laboratories, Inc. 

18. R.E.M. Hodges, D. L. Drummond, and A. Gallagher, Phys. 

Rev. A6, 1519 (1973). 


85 



19. O.S. Heavens, J. Opt. Soc. Am. 

20. J. Holstein, Phys. Rev. 72 , 1212 

21. A.J. Palmer, J. Appl. Phys. 47 , 


51, 1058 (1961). 
(1947). 

3088 (1976). 


86 



APPENDIX 


Section 1 


C XEFL - DISCHARGE PUMPED XEFL LASER 

REAL UASE , MASXE , MASH , M 1 ,M2,MI , K, K3X,K5X , K4 I ,NEO‘NE 
REAL NFP 
REAL J 
' REAL W 
REAL LENGTH 
REAP K'PLliS, KNI , KEL 
REAL L 

REAL MD1 , MD I ,MD 
EXTERNAL FUNC 

COMiVON/STOR/M, XE , NE ,WE? , C {f i ,CXI , XEFG , EM 1 ,EXi <EX-EFU,CEL, 
+MASF , M ASM, V , LENGTH , A , L , R , DT , E , FT . 

C CONSTANT DATA 

MASE=0. 5F6? E= I . 6E- 19; MAS XE= 131*1 840*. 5E6 
L=0 

C= ! OOOE -9 
D=. 1 

A=3 * f 4*(D/2)**2 
C M=HELIJM 

MA$M=4* i 840*. 5E6 
EXI=8.3; EX2=9.o » EXI=12.03 
EMI =22 . ? EM2=23.? EMI=24.6 
EXEFU=3 . 1 

CXl=.3E-i 6; CX2-.05E-1 6? CXI=.IIE-I6 . 

CX i I = 5E- 1 6 ; CX2I=5E-I6 

CM1=.0?5E-I 7 1 CM2= . 025E- ! 7 ; CMI=.08E-17 

CHII=5E-16? C«2I=5E-16 

CEL=5 . 5E-! 6 

CDIF= IE-1 5 


-LENGTH* I .6 
K= 1 4AQE I 9 
NEO= I E 1 0 
TG0=. 03 

AX2 = 4E7 ; AV2=32E / ; AXED I =4E7 ; AMD I =32E 7 ; AXEF=2E7 ; AMX= 1 0E7 * 

0 INITIAL CONDITIONS - 
WRITE (5, SO) 

50 FORMAT! 2X , 1 HM,5X-,2HFL, 5X , 2HXE, 5X , 2HV0, 5X , 1 HR , 5\ , 2HDT , 
+5X,6HIPRINT,5X„4HTMAX,5X,2HVS) 

ACCEPT*, N,FL,XE, VO, R,DT,‘IPRINT, TM AX, VS 

TG=TGO 

V=VO 

WRITE (5, 80) . 

80 FORMAT ( 4X , 4HT TME,6X, 5HVD I SC , SX , 1 H J ; 9 X, 2 HNE, 8X , 3HFN I . 

+7X,2HMI ,8X,3HXEI ) 

WRITE (5 , 85) - 

85 FORMAT ( 7X , 2HFL , 8X , 3HXE 1 , 7X , 4HXEFU , 6X, 4HXEFG, 6X, 4HGAIN,oX , 2HTG, 
+8 X , 2HTI : ) 
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C TIME STEP ; 

„ 20 ISTFP=ISTEP+T 
T-T+DT 

' • IF(T-TMAX) 21,2.1,22' 

21 - CONTINUE •' .V : . 

C C IRCU ITRESPONSE ' 

w E= N E 0+M I+X E ■! +MD T+-'X E D I + X E M I—FNT.' 

V 0= V 0 -DT * J'* A/C 
IF (VO-VS) 4,5,5 
5 V=VO; GO TO 6- -'-P • 

4 V=VS 
’6 CONTI NUF 

IERR=0? IJ'ID-1 ?EPS=0.02. 

IF ( I STEP- 1 ) . 10,10,8 
8 IF (TE— . 5 ) 9,9,10 
9 . TE= . 5 ; GO TO 2 
10 CONTINUE ’ 

. . -TEfSOLN (IND., EUNC , .,4, TO. ,EPS,'iERLU ■ 

‘ IF(IERR) 1,2,1 - 

• 1 TE=. 4 * 

.2 CONTINUE . \ - . « ' > 

C RAJE CONSTANTS; - ELCTROMIC „ •' . 

V D= S 0 RT C 2**TE /M AS M > *3 E TO ' *" ’ * 5 * . 

J=NE*E*V,D 

VD I SC= V-R*J*A-l> A*E*V D*C N E-NEP ) /DT 
NEP=NE . , •: ‘ ' * - " „ 

RC=6.38E7*CSQRT(2*TE))±*3 . v - 
Rl=RC*EXP(-EMl/TE)*t 1 .+EMl/(2*TE) ) * ( NE*CM 1 > * 

R1B=(R1/M>*EXP(EM1/TE)*M1*.'2.5 ' ' •'* \ 

R 1 X=RC*EXP( -EX 1./TE ) * (-1 . +E X 1 / ( 2 *TE ) ) * C NE*C X UXE > 

R1XB=(R1X/XFJ*EXP‘(EXI/TE)*XE'1*.25 1 ” ' ' . 

R37=RC*EXP(-EM2/TE)*C l . +EM2/C2*TE)3 *(NE*CM2*M)' V ’ 

R3 7 B= ( R3 7/M ) *EXP ( EM2/TE-) *M2* . 1 W*"' 

R37X=RC*EX.P( -EX2/TE )*( 1 . +EX2/T-2*TE) ) *(NE*CX2-OCE> 
R37XB=( R3‘7X’/XE)*EXP (EX2/TE )*XE2*. t - ' . '• " 

R3 8= RC*E XP ( t ('EM I -EM 1 > /TE) * ( t . * ( EM I -EM l ) / C 2 *TE ) ) * 

+ (NE*CM1 I'*M1 > - - - . 

R39=RC*EXP(-(EMI-EM2:)'/TE>*( 1 .+(<EMI-EM2>/(-2*TE> ) * 

+ CNE*CM2I*M2) '* > **' *-* 

R38X=RC*EXi?( -( EXI-EXt )/TE>*( 1 . +( EXI-EXH A'2*TE>) 

+* ( NE*CX 1 1 *XE 1 ) ' - . - 1 ' 

R39X=RC*EXP ( r ( EX I-EX2 )/TE ) * U +C EXI-EX2 ) /C2*TEX-> 
+*(rtE*CX2I*XE2) - ' 

R 2= RC*E XP (-EM I-/TE » (' 1 1 REM I / (2*TE> ) * ( N 6*CMT*M ) 

R 2X= RC*E XP X-EX I /TE ) * ( 1 ,+EXI/-( 2*TE ).)-*(NE*GXl'*XE) ' 
R4=,29E-8*MD.I*NE/<SQRT(TE>:> jP,':- \,1 /' 

R4X=.346E-6*XEDI*NE/(S0RT(TE)>'': ■ 1 

R4MX= i 346E-7 *XEM I *NE/,( SQRTCTE’P) ' 4? -■ 

R 1 7 = '9 . 1 E-t 0*TE**( -2.16) *EXR <*>-.' 1 65 /TE) *NE*fL 
R28=2E-8*‘XEFU*NE ‘ 

R28B=2Er-8*EXP(-EXEFU/TE ) *NE*XEFG . . 

REL= (NE*CEL* M*5G RT ( 2*TE/M ASE'>'> *3;E 10* ( 2*M ASE/M ASM >-*TF 


^^RODUGIBIblTy* OF TH@ 
O'R-GlTl AL PAGE IS POOR) 
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C RATE CONSTANTS - GAS KINETIC 
SQRTG=SORT(TG) 

K3X= ! ,E-23*EXP(. 6/TG) 

K5X= 1 ,E-23*EXP( 1 ,5/TG) 

K 4 1= I .E-23*EXP( .074/TG) 

G3«M*3. 7E-3 1 *M*MI*SQRTG 
G3X=XE* 1 2E-3 1 *M*XEI*SORTG 
G3MX=M+5 . 8E-3 ! *M*XE I +SQRTG 
G3XB= 1 ?E-3 1 *XE*XED I /K3 X+SORTG 
G5=M* ] 2E-33*M*M I *SQRTG 
G5X=M*1 5E-32*XE*XE!*S0RfG 
G5MX=M*47 E-33*M*XE I *SQRTG 
G5XB=I5E-32*XE*XED1 /K5X+SQRTG 
GI3*29E-lO*M1*XE*SQRTG 
G!4=29£-10*MD1*XE*S0RTG * 

Gi3=0; G18X=34.7E-I0*S0RTG*XE1+FL 

G20=0 ; G20X=35E- 1 O+XED 1 *FL*SORTG 

G 22= 5E-7 * M I *E N I *SO RTG 

G22 X= 1 7E-7*XE I *FN I *SORTG 

G23=3E-25*MI*FMI*M*SORTG 

G23X- I2E-25*XE I*FNI*M*SQRTG 

G24=5E-7*HUI*Fi'JI*SQRTG 

G2ax= 1 7E-7+XEDI *FN I*SORTG 

G25=3 E-2 5* ML) I*FNI *M*SORTG 

G25X= 1 2E-25*XEDI*FNI*M*SQRTG 

G27X-46E- 1 0*XEFU*FL*SOR TG 

G3I=10E-10*M1**2*SQRTG 

G3 1 X=30E- 1 0*XE 1 **2*SQRTG 

G3 I MX=20E-1 0*XE 1 *Mi*SQRTG 

G32= 1 OE- 1 0*M I *MD 1 *SORTG 

G32X=30E- ! 0*XF 1 *XED 1 *SO RTG 

G32MX«20E-10*M1*XED1*S0RTG 

G33= 1 OE— 1 0*MD ! **2*SORTG 

G33X=30E-10*XED1**2*SQRTG 

G33MX=20E-1 0*XEM 1 **2*SQRTO 

G34= 1 OE- i 0*M 1 *XEFU*SQRTG 

G34X=30E-1 0*XE 1 *XEFU*SORTG 

G34MX=0 

G35=l nE-10*MD1*XEFU*SQRTG 
G35X-30E-1Q*XEDI*XEFU*SQRTG 
G35MX=20E- ! 0*XEM I *XEFU+SORTG 
G40=F*1 E-32*F*M*SORTG 
G4l=M*2E-32+XE*F*SGRTG 
G41 B=2E-32*.M*XEFG/K41 
G44=6E-3 2*F*XE 1 ~M*SQRTO 
CSTIM=2. 6E-1 7 
C DIFFUSIN RATES 

KPLUS=. 44E3*(3E1 6/< XE+M > )*SQRT( .03/TG) 

Kit I=2E3 * (3E 1 0/ ( XE+M ) ) *SORT ( . 03/TG ) 

KEL=VD/ ( VDISC/LENGTH) 

DPLUS= 1 2* (3E 1 6/ ( XE+M ) ) *SORT(TG/.03) **3* ( 4. 3/D) **2 
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DNl=30*(3Elo/(XE+M) )*S0RT!TG/.03)**3*U.8/D)**2 
DEL=SQRT ! 2*TE/MASE ) *3E 10/(3 *CEL* ( XE+M ) ) * ( a. n/0 ) **2 
n XE=SORT ( 2 *TG/W AS XE > *3 E 1 0/ 1 3* ( M+XE ) *CD T F ) * < 4 . 8/D > *+2 
DM= DXF*SQRT ( MA SXE/M ASM ) 

DAPOS=DPLUS-KPLUS*(DPLIJS*XEI~D.iI*FNI-nEL*ME)/ 

+ ( KPI.US*XE I +KN T*FW I+KEL+NE ) 

DANEG=DN I +KN I* ( DPLUS*XE I -DN I *Fil I -DEL*ME ) / ( KPLUS *XE I+KU I*Fh T + 
+KEL*NE> 

c HATE equations 

M 1 =?M +DT*( R 1 +R4+G22+G5B-G5-G 1 8-R 1 B-313-G3 1 I ) 

M2=M2+DT*( R37+H4-AM2*M2-R3 IB) 
«I=MI+DT*(R2+G3B+R33+R39-G3-G22-O23-0A |, OS* , *I ) 

XE 1 = XF. i +DT*( R 1 X+R4*+G22 X+G5XB-G5X-G I 8X-R 1 XR-G3 i X-R3*X-JXE*XE 1 ) 

XE2= XF2 +DT* C R3 7 X+R 4X- AXE 2*XE 2- R3 7 B ) 

XEI=XEI+DT*(R2X+G3XB+G3HXB+G13+G1 4+R38 <+R39X-G3X-'33MX-G22X-G23X 
+-DAPOS+XFI) 

MDI=^DI+DT'V(G3+G31+G32+G33+G34+G3B-SM-G?A-G25-G3£i-0AFOS^^DI ) 

XED 1 = XFD 1 +DT* ( G5X+G24X-AXED* XED 1 -G5XB ) 

XEDI = XEDI+DT*(G3X+G3J X+G32X+G33X+G34X+r,3 5X-R4X»G24X-o25X-G3XB 
+-DAPOS*XEDI ) 

XEM 1 =XEM I +0T*(G5MX+G24MX-AMX*XEM I -G5MXB) 

XEMI=XEMI+DT*(G3MX+G31'. i JX+G32MX+G33MX+G34MX+G35 , 4X-R4MX-G24MX-G25MX-- 
+G3MXB-DAP0S*XEMI ) 

FNI-FNI +0T*< R 1 7-G22-G22X-G23-G23X— G24— Q24X-G25— G25X-DAwh'G*FN I ) 

FI =FL+DT* ( R 1 7B+U40-R ! 7-G 1 8-G 1 BX-G20-G20 .< ) 
XEFG«XEFG+0T*(AXEF*irXEFU+R28+G27X+G35+G41+G35X- , ?283-G41B) 
XEFU=XEFU+DT*(G1 8X+G20+G23X+G25X+R23B-AXEF*XEFU-G27X-R28-G34-G34X- 
+G35-G35X) 

F=F+DT*(RI7+G18X+G20X+G22+G22X+G24+G24X+G25+G25X+?+G27X+G34+G34X+G 
+4 IB 

+ -Q44-G4 1 -2 *G40 ) 

G <^AS TEMPERATURE 

TG=TG+DT*2/(3*M)*! HE+VD*! VD ISC/LENGTH) -( 4. 3/D) **2*K*(TG-TG0 > ) 

C GAIN 

GAIN=CSTI «*< XEFU-XEFG+EXP ( 974/TG ) ) 

C uUTPUT 

IF < ISTEP/IPRINT-FLOAT(ISTEP)/FLOAT(IPRlUT) ) 3! ,30,31 

30 CONTINUE 

WRITE (5,90) T, VDISC , J , "!E ,FN I , M I , XE I 
90 FORMAT ( 7 ( E 1 0 . 2 ) ) 

WRITE! 5, 05) FL,XEt ,XEFU,XEFG,GAIN,TG,TE 
95 FORMAT (3X,7(E10,2)) 

31 CONTINUE 
GO TO 20 

22 WRITE (5, 200) 

200 FORMAT! 2X, I9H>0, GO ON <0, END) 

ACCEPT*,/ 

IF ( Z ) 41,41,42 
42 WRITE (3, 250) 

250 FORMAT! 2X,28HNEW VALUES: DT, TMAX, I PRINT) 

A CCEPT* , DT , TMA X , I P 3 1 NT 
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ISTEP=0 
GO TO 20 
STOP 
41 ENID 

FUNCTION FUHC(TE) 

REAL M, M 1 , m2 , NE , MASE , MASM , LENGTH , NEP , L 
COMMON/STOR/M, XE , NE ,!IEP , CM ] ,CXi ,XEFG, EM I , EX1 ,EXEFU,CEL, 
+MASE , MASM , V , LENGTH , A ,L , R , DT , E , FI. 

VD=SO RT ( 2 *TE/M AS M ) *3 E 1 0 
RC=6.38E'/*(SORT(2*TE) > **3 
IF (TE-.4) 51,52,52 

51 R 1=0. ; R 1 X=0 . 5 GO TO 53 

52 CONTINUE 

R 1=RC*EXP ( -E M 1 /TE ) *( 1 . +EM 1 / ( 2*TE ) )* ( CM 1 *M ) 
R1X=RC+EXP(-EXI/fE)*( 1 .+EXI/(2*TE))*(CX1*XE) 

53 CONTINUE 

R17=9.1E-10*TE**(-2. 15>*EXP(. I65/TE)*FL 
R28B=2E-8*EXP(-EXFFU/TE)*XEFG 

REL= ( CEL*M*SQRT( 2*TE/M ASE ) ) *3E 1 0* ( 2*MA5E/MASM ) *TE 
F UNO (V-P*A*E+NE*VD-L*A*E*VD*(NE-NEP)/nT) /LENGTH 
+ -( l/VD)*(R1X*EXI+R1*E¥l+R17*TE+R28B*EXEFU+REL) 

RETURN; END 
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Section 2 


C KXLCP - PULSED K-XC LASER 

REAL M A S E f r. AS X L , M AS te , M 1 , rt 2 , M I , K , K 3 X , K 5 X f K 4 1 , U £ 0 , N E 

REAL LAMDE , L AMDC , K AE X , K AD 1 ,KXEX,KXDI 

REAL LLP 

REAL J 

KbAL LENGTH 

REAL KPLUS ,RNI ,KEL 

REAL L 

EXTERNAL FUNC 

COMilON/ STOR/XE , P0 ,CP1 ,CXl ,NE , t.’E P , L , A ,R , DT , E , 

+ LP1 ,EX1 , CEL, ilASE,MASX£,V, LENGTH , Pi ,CP1T ,CPI ,EPI 
COEttOK/STE/lL 
C CuwSTANl DATA 

HAS E= 0 . 5 E 6 ; E=1.6E-19; MASX£«131*1840* . 5E6 
R=0 

C=2 . 7E-9 
D=1 
1\-21 
C XL=XLNGR 

LAwDE= . 85E-4 
LAMDu=l . 04E-4 
AP=3 . 69E7 
RfcX=3 .7E-8 
RDI=5.4E-8 
RDI£=4 .4E-8 
DRDFE= . 53E-8/ (579) 

DRDFD=1Lt8/ ( . 04E4) 

DRDFDE= . 4E-8/ (1.388E3) 

VEA=- .079 
VEX=.G62 
VDA=- . 66 
VDX=~ .262 
VDX£=- . 56 0 
EP1=1.61 
EP2=3 . 5 
EPI=4 .34 
EPXl = 1 .46 
EPXI = 4 .19 
EPD1-1.46 
EPDI=4 .19 
CP1=6E-15 
CP2=6E— 1 6 
CPI= . 5E-16 
CP1I=5E-16 
CPXl=oE-15 
CPXI-.5E-16 
CPD1=6£-15 
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CPDI= . 5E-16 

CPHI= . 0 5E-18 

CEL=5 . 5E-16 

CDIF-1E-15 

SA=8 . 49E4 ;SB=7 .18 

LENGTH=1 

K=160E19 

AX2=4E7; AP=3 .69E7 ?AXEDl=4E7 
CSTf4E=AP/2*LAMDE* * 2*REX**2*DRDFE 
CSTMD-AP/ 2* LAMDD* * 2*RDI ** 2*DRDFD 
CSTMDE=AP/2*LAMDE**2*RDIS**2*DRDFDE 


C INITIAL CONDI'XIOnS 
aKITE ( 5 , 50) 

50 FORMAT ( 2X f 2HXE , 2X , 2RTG , 2X > 1HW , 2X , 1HV/2X , 1HC , 2X , 2EDT , 2X , 6HIPRINT , 
+2X , 4H1MAX , 2X , 5HDELAY) 

ACCEPT* , XE , TG , w , V ,C , DT , IPRINT , TMAX , DELAY 

Pfr=2 . 7E16* (27 3/(27 3+TG ) ) *10** ( SE- . 05 2*SA/ (TG+27 3) ) 

TG= (TG+27 3) * (1/1.16E4) 

TG0=TG 

F=W/ (1,6E-19*EPI)*(1. 24/SPI- . 2) 

KXDI=1 . 8E-22*EXP ( ,56/TG) 

KXEX=2 . 2E-23*EXP (— . 05/TG) 

PX=KXEX*P0*XE 
PD=KXDI*P0**2 
PX=KXEX*P0*XE 
-- ' PI=F*GPHI *P3* DELAY - 
Pi— 10 . *PI 

AT=AP*1.6/(8.3E-15*P0*D) 

- - IF'fAT— AP) "18,17,17- - 

17 AT=AP 

18 CONTINUE 

- - WRITER-, 82)- P0,PD,PI 

82 FORMAT (2X, "P0= ' , El 0 . 3 , 5X , 'PD= ' , E10 . 3 , 5X , 'PI**',E10 .3) 

WRITE (5, 80) 

8 0 - -FORMAT ( 4X-,-4 RTI ME , 6X , 5HVD7SC , 5X , 1HJ , 9X-, 2HNE , 8X , 2 HP I , - 

+8X,3HPX1,7X,3HPD1) 

WRITE (5, 85) 

— - 8 5 FORMAT ( 7X , 1HP , 9X , 2HP1 , 8X ,-5HG AINE , 5X ,-5HGAIND , 5X , 5HBETAE , 5X , 5HBETAD , 

+5X, 2HTE) 

C TIME STEP 

20 ISTEP=ISTEP+1 
T=T+DT 

IF (T-TMAX) 21,21,22 

21 CONTINUE 

C CIRCUIT RESPONSE 
P=P0- (Pl+PI) 

NE=PI 

V=V-DT* J* A/C 

IERR=t) ) IND=1; EPS=0 .01 

TE=SOLN ( IND ,-FONC ,.1,5. ,EPS,IERR) 
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IF(IERR) 1,2,1 
1 WRITE (5,*) IERR;GO TO 30 

STOP 

.2 CONTINUE 

C KATE CONSTANTS - ELCTRONIC 

VD=SQRT(2*TE/MASXE) *3E10 
J=NE*E*VD 

VDISC=V-R* J*A-L X A*E* VD* (NE-NfiP) /DT 
N£P=NE 

RC-6 . 38£7* ( SQRT ( 2*TE) ) **3 

RP1=RC*EXP (-EP1/TE) * (1+EPl/ (2*TE) ) * (NE*CPl) *P 

RPlB=RPl/P*EXP (EP1/TE) *Pl*.33 

RPI=RFN (EPI) * (NE*CPI) *P 

RP1I=KFN (EPI-EP1) * (NE*CPlI ) *P1 

RPXl=RFN ( EPX1 ) * (NE*CPXl) *PX 

RPX1B=KPX1/PX*EXP (EPX1/TE) *PX1* .5 

RPXI=KFN (EPXI) * (N£*CPX1) *PX 

RPDl-RFN (EPD1) * (NE*CPD1) *PD 

RPDlB=RPDl/PD*EXP (EPDl/TE) *PDl*l . 

RPDI=RFN (EPDI ) * (NE*CPDI) *PD 
RLCP=PI*2 . 23E-26*NE*TE** (-4 .39) *NE 
RECPX=3 . 3E-5*TE** <- . 67) *PXI*NE 
RECPu=3 . 3E~5*T£** . 67) *PDI*NE 
C RATE CONSTANTS - GAS KINETIC 
SQRTG=SQRT(TG) 

KAEX=2 . 2 £-2 3* EXP ( .074/TG) 

KADI=6.5E-23*EXP( .735/TC-) 

KXEX-2 . 2E— 23*EXP {- . 05/TG) 

KXDI=1 . 6E-22*EXP ( . 56/TG) 

GPX=XE* 8E-32*XE* Pi 
GPXB=XE*8E— 32/KAEX 
GPD=P* (XB*8E-3E*Pl) 

GPDB= (XL*8E-30) /KADI 
GPxI= (X£*8t-32*XE) *PI 
GPX1B= (XE*8E-32*PXI ) / KAEX 
GPDI= (XE*8E-3w*P) *PI 
GPDI B= ( XE*8E~3fe5*PDI ) /KADI 
C DIFFUSION KATES 

KPLUS = 2E3* (3E16/(XE) ) *SQFT ( . 83/i’G) 

KbL=VD/ ( VDISC/LENGTB) 

DXE=SQRT ( 2*TG/fiASXE) *3E10/(3* (XE) *CDIF) * (4 .8/D) **2 
DAPOS=KPLUS*TE* (4.8/D) **2 
DP=DXE*1 . 8 
C RATE EQUATIONS 

-- - - Pl=Pl+DT* (KP1+RECP+RECPX+RECPD-KP1B-AP* (PXl+PDl) -DP*P1) 

IF (Pl/P-3*EXP (-EP1/TE) ) 8,9,9 
9 P1=P*3*EXP (-LP1/TE) 

- 8 CONTINUE 

PX=KXEX*P*XE 

PD=KXDI*P**2 

-PI=PI+DT* (RPI+RPl-I+RP2l+GPN*P+P*CPHI*F-RECP-RECPX-RECPD-DAPOS*PI) 
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PXl=GPX/ (GPXB+AP) 

PXI=KAEX*PI*XE 

-PBl-=GPD/-fGP D B - f AP) - - 

PDI=KADI*P*PI 
C GAS TEMPERATURE 

TG=4FG+DT*-{NE*VD*-VDISe/LENGTH-{4 .8/D) **2*K* (TG-TG0)') *'2/{3*XE) 

C GAIN 

B£TAE=CSTME*2*P*XE*EXP(-VEX/TG)/3E10+CSTMDE*P**2*EXP(-VDXE/ 

KEG-) 25/3E10 - - - - - 

BETAD=CST-MD*P**2*EXP (-VDX/TG) * . 25/3E10 

GAINE=CSTME* ( PXl*EXP (-VEA/TG) / (KAEX*1 .5) -2*P*XE*EXP (-VEX/TG) ) /3E10 

+-eST?MBE*-{P**2*E'X'P-(-VDXE/TG)*.2 5)/-3El0 -- - — 

GAIND=CSTMD* (PDl*EXP (-VDA/TG) / (KADI *12) -P**2* EXP (-VDX/TG) * . 25 ) 
+/3E10 

- C -OUTPUT - - — 

IF (ISTEP/IPRINT-FLOAT ( ISTEP) /FLOAT { IPRINT) ) 31,30,31 

30 CONTI MG E 

WR-IT-E ( 5 , * ) TG ... 

WRITE{5 ,90) T,VDISC,J ,NE,PI ,PXl ,PDl 
90 FORMAT (7 (E10.2) ) 

. WRITE ( 5 , 95) - P ,Pi fGAINE ,GAIND , BETAS , BETAD ,TE 

95 FORMAT ( 3X , 7 (E10 . 2) ) 

31 CONTINUE 
- -GO— TO *20 

22 WRITE (5, 203) 

200 FORMAT (2X,19H>0 , GO ON <0 , END) 

• - -ACCEPT* , 2 

IF (2 ) 41,41,42 
42 WRITE (5 , 250) 

250 FORMAT (2X,28HNEW VALUES: DT> TMA-X , IPRINT) 

ACCEPT* , DT ,TMAX , IPRINT 
ISTEP=B 
GO TO 20 
STOP 
41 END 

FUNCTION FUNC(TE) - - -■ 

REAL ME, MASE,MASXE, LENGTH, NEP,L 
COMMON/STOR/XE , P 0 , CPI , CXI , NE , NE P , L , A , R , DT , E , 

+ EPl ,EX1 ,CEL ,MASE ,MASX£ ,V , LENGTH ,Pl ,CPlI ,CPI ,EPI 
VD=SQRT (2*TE/MASXE) *3E10 
RC=6 . 38E7* (SQRT ( 2*TE) ) **3 
RPl=RC*EXP (-EP1/TE) * (1+EPl/ (2*TE) ) *CP1*P0 
REL=(CEL*XE*SQRT(2*TE/MASE) ) *3E10* { 2*MASE/MASXE) *TE 
FUNC= (V-R*A*E*NE*VD-L*A*E*VD* (NE-NEP)/DT) /LENGTH 
+ — (1/VD) * (RPl*EPl+REL) 

RETURN ; END 
FUNCTION RFN(Y) 

COtoMON/STE/TE 

RFN=6 .3807* (SQRT (2*TE) ) **3*EXP (-Y/TE) * (1+Y/ (2*TE) ) 

RETURN; END 


EX$ 

OUTPUT FILL: KXETV.FOR;15 [New version] 



